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Abstract: Jet shapes are weighted sums over the four- momenta of the constituents of a 
jet and reveal details of its internal structure, potentially allowing discrimination of its par- 
tonic origin. In this work we make predictions for quark and gluon jet shape distributions 
in iV-jet final states in e + e~ collisions, defined with a cone or recombination algorithm, 
where we measure some jet shape observable on a subset of these jets. Using the framework 
of Soft-Collinear Effective Theory, we prove a factorization theorem for jet shape distri- 
butions and demonstrate the consistent renormalization-group running of the functions in 
the factorization theorem for any number of measured and unmeasured jets, any number 
of quark and gluon jets, and any angular size R of the jets, as long as R is much smaller 
than the angular separation between jets. We calculate the jet and soft functions for angu- 
larity jet shapes r a to one-loop order [0{a s )) and resum a subset of the large logarithms of 
T a needed for next-to-leading logarithmic (NLL) accuracy for both cone and ki-type jets. 
We compare our predictions for the resummed r a distribution of a quark or a gluon jet 
produced in a 3-jet final state in e + e~ annihilation to the output of a Monte Carlo event 
generator and find that the dependence on a and R is very similar. 
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1. Introduction 

1.1 Motivation and Objectives 

Jets provide troves of information about physics within and beyond the Standard Model of 
particle physics. On the one hand, jets display the behavior of Quantum Chromodynamics 
(QCD) over a wide range of energy scales, from the energy of the hard scattering, through 
intermediate scales of branching and showering, to the lowest scale of hadronization. On 
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the other hand, jets contain signatures of exotic physics when produced by the decays of 
heavy, strongly-interacting particles such as top quarks or particles beyond the Standard 
Model. 

Recently, several groups have explored strategies to probe jet substructure to distin- 
guish jets produced by light partons in QCD from those produced by heavier particles 
[1, 2, 3, 4, 5, 6, 7, 8], and methods to "clean" jets of soft radiation to more easily iden- 
tify their origin, such as "filtering" or "pruning" for jets from heavy particles [5, 9, 10] or 
"trimming" for jets from light partons [11]. Another type of strategy, explored in [12], to 
probe jet substructure is the use of jet shapes, which are modifications of event shapes [13] 
such as thrust. Jet shapes are continuous variables constructed by taking a weighted sum 
over the four-momenta of all particles constituting a jet. Different choices of weighting 
functions produce different jet shapes, and can be designed to probe regions closer to or 
further from the jet axis with greater sensitivity. 1 While such jet shapes may integrate over 
some of the detailed substructure for which some other methods search, they are better 
suited to analytical calculation and understanding from the underlying theory of QCD. 

In this paper, we consider measuring the shape of one or more jets in an e + e~ collision 
at center-of-mass energy Q producing N jets with an angular size R according to a cone or 
recombination jet algorithm, with an energy cut A on the radiation allowed outside of jets. 
We use this exclusive characterization of an iV-jet final state looking forward to extension 
of our results to a hadron collider environment, where such a final state definition is more 
typical. For the jet shape observable we choose the angularity r a of a jet, defined by (cf. 



where a is a parameter taking values — oo < a < 2 (for IR safety, although factorizability 
will require a < 1), the sum is over all particles in the jet, Ej is the jet energy, 
is the transverse momentum relative to the jet direction, and r/ = — lntan(#/2) is the 
(pseudo)rapidity measured from the jet direction. The jet is defined by a jet algorithm, such 
as a cone algorithm, the details of which we will discuss below. We complete the calculation 
for the jet shape r a for jets defined by cone or recombination algorithms, but our logic and 
methods could be applied to a wider spectrum of jet shapes and jet algorithms. We have 
organized our results in such a way that the pieces independent of the choice of jet shape 
and dependent only on the jet algorithm are easily identifiable, requiring recalculation only 
of the observable-dependent pieces to extend our results to other choices of jet shapes. 

Reliable theoretical prediction of jet observables in the presence of jet algorithms is 
made challenging by the presence of many scales. Logarithms of ratios of these scales can 
become large and spoil the behavior of perturbative expansions predicting these quantities. 
These scales are determined by the jet energy u), the cut on the angular size of a jet R, 
the measured value of the jet shape such as r a , and any other cut or selection parameters 
introduced by the jet algorithm. 

1 The original "jet shape," to which the name properly belongs, is the quantity ^(r/R), the fraction of 
the total energy of a jet of radius R that is contained in a subjet of radius r [14, 15, 16] . This observable 
falls into the larger class of jet shapes we have described here and for which we have hijacked the name. 



[12, 17]) 




(1.1) 
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Precisely this separation of scales, however, allows us to take advantage of the powerful 
tools of factorization and effective field theory. Factorization separates the calculation of 
a hard scattering cross section into hard, jet and soft functions each depending only on 
physics at a single scale [18, 19]. Renormalization group (RG) evolution of these functions 
between scales resums logarithms of these scales to all orders in a s , with the logarithmic 
accuracy determined by the order to which the anomalous dimensions in the running are 
calculated [20] . Effective field theory organizes these concepts and tools into a conceptually 
simple framework unifying many ingredients going into traditional methods, such as power 
counting, gauge invariance, and resummation through RG evolution. The rules of effective 
theory facilitate proofs of factorization and achievement of logarithmic resummation at 
leading order in the power counting and make straightforward the improvement of results 
order-by-order in power counting and logarithmic accuracy of resummation. 

1.2 Soft-Collinear Effective Theory and Factorization 

Soft-Collinear Effective Theory (SCET) [21, 22, 23, 24] has been successfully applied to the 
analysis of many hard scattering cross sections [25] including the production of jets. SCET 
is constructed by integrating out of QCD all degrees of freedom except those collinear to 
a lightlike direction n and those which are soft, that is, have much lower energy than the 
energy of the hard scattering or of the jets. Using this formalism, the factorization and 
calculation of two-jet cross sections and event shape distributions in SCET were developed 
in [26, 27, 28, 29]. Later, these techniques were extended to the factorization of jet cross 
sections and observables using jet algorithms in [30]. Calculations in SCET of two-jet rates 
using jet algorithms have been performed in [27, 31], and more recently in [32]. Calculations 
of cross sections with more than two jet directions have been given in [33, 34, 35]. 

Building on many of the ideas in these previous studies, in this paper, we will demon- 
strate a factorization theorem for jet shape distributions in e + e~ —> N jet events, 

^(Pi,.. Pjv) = ff(0) p N)H ( nitUl . . . . nN ^ N . M) (L2) 

OT\ ■ ■ ■ UTm 

X Jm,uii(n; M) " " " Jn M ,ui M (TM; H) <g> S ni ... nN (Ti, . . . ,TM',R, A;//) 
X ^™M+i,^A/+i \R'i A 4 ) ' ' ' Jn N ,u) N {R\ A*) j 

where the TV" jets have three-momenta Pj, and M < N of the jets' shapes t±, . . . ,tm are 
measured, is the Born cross-section, H is a hard function dependent on the directions 
rii and energies Wj of the N jets, J n ^ (r) is the jet function for a jet whose shape is measured 
to be r, J n ,u)(R) is the jet function for a jet with size R whose shape is not measured, and 
S is the soft function connecting all N jets, dependent on all jets' shapes r», sizes R, and 
total energy A that is left outside of all jets. The symbol "(8)" stands for a set of convolution 
integrals in the variables Tj between the measured jet functions and the soft function. All 
terms in the factorization theorem depend on the factorization scale /i. 

SCET is typically constructed as a power expansion in a small parameter A formed by 
the ratio of soft to collinear or collinear to hard scales, determined by the kinematics of the 
process under study. A is roughly the typical transverse momentum pt of the constituent 
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of a jet (relative to the jet direction) divided by the jet energy Ej. This is set either by 
the measured value of the jet shape r a for a measured jet or the algorithm measure R 
for an unmeasured jet. Thus we encounter in this work the new twist that the size of 
A may be different for different jets. We will comment on further implications of this in 
subsequent sections. Still, in each separate collinear sector, the momentum p n of collinear 
modes in the light-cone direction n in SCET is separated into a large "label" momentum p n 
containing O(Ej) and O(XEj) components and a "residual" component of 0(X 2 Ej), the 
same size as soft momenta. Effective theory fields have dynamical momenta only of this 
soft or residual scale. This fact, along with the fact that soft quarks and soft gluons can be 
shown to decouple from collinear modes at the level of the Lagrangian [24] , makes possible 
the factorization of a jet shape distribution into hard, jet, and soft functions depending 
only on the dynamics at those respective scales. 

In using SCET for jets in multiple directions and using jet algorithms to define the 
jets, we will encounter the need for several additional criteria to ensure the validity of the 
iV-jet factorization theorem. 

• First, to ensure that the algorithm does not group final-state particles into fewer than 
N jets, the jets must be "well separated." This allows us to use as the effective theory 
Lagrangian a sum of N copies of the collinear part of the SCET Lagrangian for a 
single direction n and a soft part, and to construct a basis of iV-jet operators built 
from fields from each of these sectors to produce the final state. Our calculations will 
reveal the precise quantitative condition that jets must satisfy to be "well separated" . 

• Second, to ensure that the jet algorithm does not find more than N jets, we place an 
energy cut A on the total energy outside of the observed jets. We will take this energy 
A to scale as a soft momentum so that we will be able to identify the total energy of 
each jet with the "label" momentum on the SCET collinear jet field producing the 
jet. Corrections to this identification are subleading in the SCET power counting. 

• Third (and related to the above two), we will assume that the iV-jet restriction on the 
final state can itself be factorized into a product of N 1-jet restrictions, one in each 
collinear sector, and a 0-jet restriction in the soft sector. We represent the energetic 
particles in the ith jet by collinear fields in the SCET Lagrangian in the rjj collinear 
sector and soft particles everywhere with fields in the soft part of the Lagrangian. 
We then stipulate that the jet algorithm acting on states in the rii collinear sector 
find exactly one jet in that sector, and when acting on the soft final state find no 
additional jet in that sector. 

• Fourth, the way in which a jet algorithm combines particles in the process of finding 
a jet must respect the order of steps envisioned by factorization. In particular, fac- 
torization requires that the jet directions and energies be determined by the collinear 
particles alone, so that the soft function knows only about the directions and colors of 
the jets, not the details of any collinear recombinations. Ideally, all energetic collinear 
particles should be recombined first, with soft particles within a radius R of the jet 
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axis being recombined into the jet only afterwards. Jet algorithms in use at experi- 
ments do not have this precise behavior, but we will discuss in Sec. 3.4 the extent to 
which common algorithms meet this requirement and estimate the size of the power 
corrections due to their failure to do so. In general, we will find that for sufficiently 
large R, infrared-safe cone algorithms and kx-type recombination algorithms satisfy 
the requirements of factorizability, with anti-kx allowing smaller values of R than kp. 

After enforcing the above requirements, a key test of the consistency of Eq. (1.2) will 
be the independence of the physical cross section on the factorization scale \i. This requires 
the anomalous dimensions of the hard, jet, and soft functions to sum to zero, 

o = ft* (m) + 7Jm+iG r ; m) + • ■ • + jj n (R; vMn) ■ ■ ■ 5{r M ) 

+ Ui ( T ii a0<5( t 2) ■ ' ■ <Ktzv) H h <5(ri) • ■ ■ 8(tm-i)ij m {t m ; A«) (1-3) 

+ 7s(n,..-, tat, #;//). 

It seems highly nontrivial that this condition would be satisfied for any number, size, and 
flavors of jets (and that the soft anomalous dimension be independent of A), but we will 
demonstrate that it does hold at 0(a s ), up to corrections of 0(l/t 2 ) which violate Eq. (1.3), 
where t is a measure of the separation between jets. In particular, for a pair of jets, i, j, 
with 3- vector directions separated by a polar angle ipij, the separation tij is given by 

tan(^ fc>{ /2) 

**' ,= tan^/2) • (L4) 

Now define t (no indices) as the minimum of tij over all jet pairs. This quantifies the 
qualitative condition of jets being well-separated, f > 1, that is required to justify the 
factorization theorem Eq. (1.2). The factorization theorem is valid up to corrections of 
0(A) in the SCET power expansion parameter and corrections of 0(l/t 2 ) in the separa- 
tion parameter. As an example of the magnitude of t, for three jets in a Mercedes-Benz 
configuration (ip = 2ir/3 for all pairs of jets), 1/t 2 = 0.04 for R = 0.7 and 1/t 2 = 0.1 for 
R = 1, so these corrections are indeed small. More generally, for non-overlapping jets, 
V> > 2R, we have 1/i 2 < 1 /4. 

Notice that for back-to-back jets {ip = n), t — > oo. Thus, for all cases previously con- 
sidered in the literature, the jets are infinitely separated according to this measure, and no 
additional criterion regarding jet separation is required for consistency of the factorization 
and running. A key insight of our work is that for an iV-jet cross-section described by 
Eq. (1.2), the factorization theorem receives corrections not only in the usual SCET power 
counting parameter A, but also corrections due to jet separation beginning at 0(l/t 2 ). 

1.3 Power Corrections to Factorized Jet Shape Distributions 

As always, there are power corrections to the factorization theorem which we must ensure 
are small. One class of power corrections arises from approximating the jet axis of the 
measured jet with the collinear direction 7ij, which labels that jet in the SCET Lagrangian. 
This direction rij is the direction of the parent parton initiating the jet. The jet observable 
must be such that the difference between the parent parton direction and the jet axis 
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identified by the algorithm makes a subleading correction to the calculated value of the jet 
observable. In the context of angularity event shapes, such corrections were estimated in 
[17, 29] and found to be negligible for a < 1, and we find the same condition for jet shapes. 

In the presence of algorithms, however, there are additional power corrections due to 
the difference in the soft particles that are included or excluded in a jet by the actual 
algorithm and in its approximated form in the factorization theorem. We study the effect 
of this difference on the measurement of jet shapes, and find that for sufficiently large R the 
power corrections due to the action of the algorithm on soft particles remain small enough 
not to spoil the factorization for infrared-safe cone and kp-type recombination algorithms. 
Algorithm-related power corrections to jet momenta were studied more quantitatively in 
[36], and their estimated R dependence is consistent with our observations. 

We do not address in this work the issue of power corrections to jet shapes due to 
hadronization. Event shape distributions are known to receive power corrections of the 
order l/(r a Q), enhanced in the endpoint region but suppressed by large energy. The 
endpoints of our jet shape distribution near r a — > 0, therefore, will have to be corrected by 
a nonperturbative shape function. Such functions have been constructed for event shapes 
in [37, 38]. The shift in the first moment of event shape distributions induced by these 
shape functions was postulated to take a universal form in [39, 40] based on the behavior 
of single soft gluon emission, and the universality was proven to all orders in soft gluon 
emission at leading order in the SCET power counting in [41, 42]. This universality relied 
on the boost invariance of the soft function describing soft gluon radiation from two back- 
to-back collinear jets. The extent to which such universality may survive for jet shapes 
with multiple jets in arbitrary directions is an open question that must be addressed in 
order to construct appropriate soft shape function models to deal adequately with the 
power corrections to jet shapes from hadronization. Nonperturbative power corrections to 
jet observables from hadronization and the underlying event in hadron collisions were also 
studied in [36], and hadronization corrections were found to scale like 1/R. In this work, 
we focus only on the perturbative calculation and resummation of large logarithms of jet 
shapes, and leave inclusion of nonperturbative power corrections for future work. 

1.4 Resummation and Logarithmic Accuracy 

Knowing the anomalous dimensions of the hard, jet, and soft functions in the factorization 
theorem allows us to resum logarithms of ratios of the hard, jet, and soft scales. We 
take this opportunity to explain the order of accuracy to which we are able to resum 
these logarithms. For an event shape distribution da/dr (i.e. Eq. (1.2) with two jets and 
integrated against 8{t — t\ — T2)), the accuracy of logarithmic resummation [43] is typically 
characterized by counting logs in the exponent In R(t) of the "radiator," 



where they appear in the form a™ln m r with m < n + 1. At leading-logarithmic (LL) 
accuracy all the terms with m = n + 1 are summed; next-to-leading-logarithmic (NLL) 
accuracy sums also the m = n terms, and so on. In more traditional methods in QCD, event 
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shapes that have been resummed include NLL resummation of thrust in [43, 44], jet masses 
in [43, 45, 46, 47], jet broadening in [48, 49], the C-parameter in [50], and angularities in 
[17]. Resummation of an event shape distribution using the modern SCET method was 
first illustrated with the thrust distribution to LL accuracy in [51]. Heavy quark jet mass 
distributions were resummed in SCET to NLL accuracy as part of a proposed method to 
extract the top quark mass in [52] . The N 3 LL resummed thrust distribution in SCET was 
compared to LEP data to extract a value for the strong coupling a s to high precision in 
[53]. Angularities were resummed to NLL accuracy in SCET in [54] directly in r a -space 
instead of in moment space as in [17]. 

Summation of logarithms in effective field theory is achieved by RG evolution. In the 
factorized radiator of the thrust distribution Eq. (1.5), one finds that the hard function 
contains logarithms of n/Q, the jet functions contain logarithms of /j,/(Qy/r), and the soft 
function contains logarithms of fi/(Qr). Thus, evaluating these functions respectively at 
the hard scale hh = Q, jet scale \xj = Q\fr and soft scale /is = Qt eliminates large 
logarithms in each function. They can then be RG-evolved to the common factorization 
scale [i after calculating their anomalous dimensions. The solutions of the RG evolution 
equations are of the form that logarithms of r are resummed to all orders in a s to a 
logarithmic accuracy determined by the order in a s to which the anomalous dimensions 
and hard/jet/soft functions are known. This underlying hierarchy of scales is illustrated 
Fig. 1 [in this case, with only one (measured) jet scale and soft scale and ui = Q] along 
with a table that lists the order in a s to which various quantities must be known in order 
to achieve a given N fc LL accuracy in the exponent of the radiator Eq. (1.5). The power of 
the EFT framework is to organize of the logs of r arising in Eq. (1.5) into those that arise 
from ratios of the jet to the hard scale and those that arise from ratios of the soft to the 
hard scale, which then allows RG evolution to resum them. 

For the multijet shape distribution in Eq. (1.2), the strategy to sum logs is the same, 
but is complicated by the presence of additional scales. This also makes trickier the clas- 
sification of logarithmic accuracy into the standard N fc LL scheme. Our aim will be to 
sum as many logs of the jet shapes r a as possible, while not worrying about any others. 
For instance, phase space cuts induce logs of R and A/uj (where a; is a typical hard jet 
energy), and the presence of multiple jets induces logs of jet separations rij • rij or ratios of 
jet energies uji/ujj. We will not aim to sum these types logs systematically in this paper, 
only those of r a (though we sum subsets of the others incidentally). In particular, resum- 
ming the phase space logs of R or A/uj is complicated by how the phase space cuts act 
order-by-order in perturbation theory 2 , and the fact that a simple angular cut R is less 
restrictive than a small jet mass or angularity on how collimated a jet must be. That is, 
an angular cut allows particles in a jet to be anywhere within an angle R of the jet axis 
regardless of their energy, while a small jet mass or angularity forces harder particles to be 
closer to the jet axis. The former allows hard particles to lie along the edges of a jet, and 

2 The JADE algorithm is one well-known example in which resummability even of leading logarithms of 
the jet mass cut y is spoiled by the differences in the jet phase space at different orders in perturbation 
theory [55]. Another example that will not work is using a kT-type algorithm with R randomly chosen for 
each recombination. This is clearly such that resummation of logarithms of R cannot be achieved. 
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Figure 1: An illustration of generic scales along with a table of log-accuracy versus perturbative 
order. A cross section with jets of energy ~ w, radius R, and energy A outside the jets, with 
some jets' shapes r a being measured and others' shapes left unmeasured, induces measured and 
unmeasured jet scales at /Zj? eas and ^ t y imeas . Dynamics at these scales are described by separate 
collinear modes in SCET. Soft dynamics occur at several soft scales, /za and induced by the 
soft out-of-jet energy cut A and jet radius R, and /^™ eas induced by the measured jet shape r a . RG 
evolution in SCET resums logs of ratios of jet scales to the hard scale fin individually, and logs 
of the ratio of a "common" soft scale /is to the hard scale. Remaining logs of ratios of soft scales 
to one another are not resummed in current formulations of SCET. The accuracy of logarithmic 
resummation of these ratios of scales is determined by the order to which anomalous dimensions and 
matching coefficients or matrix elements (i.e. hard/jet/soft functions) are calculated in perturbation 
theory. In this paper we perform the RG evolution indicated by the arrows to NLL accuracy. 

soft radiation from such configurations that escapes the jets can lead to logs of A/uj that 
are not captured in our treatment. These are not issues we solve in this paper, in which 
we focus on resumming logs of jet shapes r a . (Some exploration of phase space logarithms 
in SCET was carried out in [31, 32].) 

A way to understand how we sum logs and which ones we capture is presented in Fig. 1. 
The factorization theorem Eq. (1.2) organizes logs in the multijet cross section into those 
in the hard function, those in measured jet functions, those in unmeasured jet functions, 
and those in the soft function, much like for the simple thrust distribution. The difference 
is the presence now of multiple jet and soft scales. Logarithms in jet functions can still be 
minimized by choices of individual jet scales, for a jet whose shape r a 

is measured, and ^^ nmeas ~ cjtan(i?/2) for a jet whose shape is not measured but has a 
radius R. Thus logs arising from ratios of these scales to the hard scale can be summed 



- 9 - 



completely to a desired N fc LL order. The complication is in the soft function. The soft 
function is sensitive to soft radiation into measured and unmeasured jets and outside of all 
jets. As we will see by explicit calculation, this induces logs of fj,tan 1 " a (R/2) / (ujT a ) from 
radiation into measured jets, and logs of fx/(2A) and /x/(2Atan^) from radiation from 
unmeasured jets cut off by the energy A. In addition, though not illustrated in Fig. 1, 
there can be logs of multiple jet shapes to one another, t^/tI. No single choice of a soft 
scale us will minimize all of these logs. 

In the present work, we will start with the simple strategy of choosing a single soft 
scale fis ~ ur a / tan 1_a (i?/2) for a jet whose shape r a we are measuring and logs of which 
we aim to resum. We will calculate hard/jet/soft functions and anomalous dimensions 
corresponding to "NLL" accuracy listed in Fig. 1. By this we do not mean all potentially 
large logs in Eq. (1.2) are resummed to NLL, but only those logs of ratios of a jet scale to 
the hard scale or of the (common) soft scale to the hard scale. We do not attempt to sum 
logs of ratios of soft scales to one another completely to NLL accuracy (which can contain 
r a ). In the case that all jets' shapes are measured and are similar to one another, r* ~ r^, 
our resummation of large logs of r* would be complete to NLL accuracy. 

We will nevertheless venture to propose a framework to "refactorize" the soft function 
into further pieces dependent on only a single soft scale at a time and perform additional 
RG running between these scales to resum the additional logarithms, and will implement 
it at the level of the 0{a s ) soft functions we calculate. However, one cannot really address 
mixed logarithms such as log(r*/ra) that arise for multiple jets until 0(a 2 ), the first order 
at which two soft gluons can probe two different physical regions. This lies beyond the 
scope of the present work. (We note, however, that our implementation of refactorization 
using the one-loop soft function does already seem to tame logarithmic dependence on A 
in our numerical studies of jet shape distributions.) 

These issues are related to some types of "non-global" logarithms described by [46, 56, 
57, 58] that spoil the simple characterization of NLL accuracy. In [59] these were identified 
as next-to-leading logs of AR 2 / (cjjT a ) and K/Q (when R <^1) that appear at 0(a 2 ) in jet 
shape distributions. These authors organized the radiator for a single jet shape distribution 
into a "global" and "non- global" part [58, 59], 

R(ri,R,A-u Ji ,Q) = R gl Ui,R,-)R ng [—,-\ . (1.6) 

In this language, the calculations we undertake in this paper resum logs in the global part 
to NLL accuracy but not in the non- global part. The first argument in R ng is related 
to ratios of soft scales illustrated in Fig. 1, and the second argument arises when there 
are unmeasured jets. In the case that all jets are measured, R ~ 1, and A ~ WjT*, the 
non-global logs vanish. 

While summing all global and non-global logs to at least NLL accuracy will be impor- 
tant for precision jet phenomenology, what we contribute in this paper are key developments 
and calculations necessary to resum even global logs of jet shapes for jets defined with al- 
gorithms. We also believe the effective theory approach and the idea of refactorizing the 
soft function will help us understand and resum many types of non- global logarithms. 
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1.5 Detailed Outline of This Work 

In this paper, we will formulate and prove a factorization theorem for distributions in the 
jet shape variables we introduced above, calculate the jet and soft functions appearing in 
the factorization theorem to 0{a s ) in SCET, and use the renormalization group evolution 
of these functions to sum global logs of r a to NLL accuracy. We consider N jets (defined 
with a cone or bp algorithm) produced in an e + e _ collision, with M of the jets' shapes 
(angularities) being measured. The key formal result is our demonstration of Eq. (1.3), 
the consistency of the anomalous dimensions of hard, jet and soft functions to 0(a s ) for 
any number of total jets, any numbers of quark and gluon jets, any number of these jets 
whose shapes are measured, and any value of the distance measure R in cone or kx- 
type algorithms (as long as t S> 1). These results lead to accurate predictions for the 
shape of the r a distribution near the peak, but not necessarily the endpoints for very 
small r a (where hadronization corrections dominate) and very large r a (where fixed-order 
NLO QCD corrections take over, which are not yet calculated and not captured by NLL 
resummation). 3 

In Sec. 2 we describe in detail the jet shapes and jet algorithms that we use. We describe 
features of an "ideal" jet algorithm that would respect exactly the order of operations 
envisioned in the factorization theorem derived in SCET, and the extent to which cone 
and recombination algorithms actually in use resemble this idealization. 

In Sec. 3, using the tools of SCET, we will derive in detail a factorization theorem 
for exclusive 3-jet production where we measure the angularity jet shape of one of the 
jets, and then perform the straightforward extension to iV-jet production with M < N 
measured jets. We will give a review of the necessary technical details of SCET in Sec. 3.1. 
In the process of justifying the factorization theorem, we identify the new requirements 
listed above on iV-jet final states and jet algorithms that must be satisfied for factorization 
to hold. In Sec. 3.4 we explore in some detail the power corrections to the factorization 
theorem due to soft radiation and the action of jet algorithms that cause tension with these 
requirements, and argue that for sufficiently large R in infrared-safe cone and recombination 
algorithms, these corrections are sufficiently small. 

Next we calculate to 0(a s ) the jet and soft functions corresponding to N cone or 
kp-type jets, with M jets' shapes measured. 

In Sec. 4 we calculate the jet functions for measured quark jets, J2(T a ), unmeasured 
quark jets, J<2, measured gluon jets, J£(r a ), and unmeasured gluon jets, J2, where ui = 2Ej 
is the label momentum of the collinear jet field in each jet function. We find that in collinear 
sectors for measured jets, the collinear scale (and thus the SCET power counting parameter 
in that sector Aj) is given by LOiTa^ 2 a \ and in unmeasured jet sectors, Aj ~ tan(.R/2). 
In studying power corrections, however, as mentioned above, we find that R must be 
parametrically larger than r a . So, in collinear sectors for measured jets, Aj is set by the 

3 Jet shapes were also studied in the QCD factorization approach in [60]. In that work QCD jet functions 
for quark and gluon jets defined with an algorithm and whose jet masses m 2 are measured were calculated 
to C(a s ). The jet mass 2 corresponds to r a for a = 0, To = m 2 /ui 2 (to < 1). A fixed-order QCD jet function 
as defined in [60] is given by the convolution of our fixed-order SCET jet function and soft function for a 
measured jet away from r a — 0. 
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shape r a with R ~ A°, while in unmeasured jet sectors, Aj ~ tan(i?/2). Thus one should 
understand tan(i?/2) to be significantly less than 1 but much larger than any jet shape r a . 

In Sec. 5 we calculate the soft function. To do this, we split the soft function into several 
contributions from different parts of phase space in order to facilitate the calculation and 
elucidate its intuitive structure. We find it most convenient to split the soft function into 
an observable-independent part that arises from soft emission out of the jets, 5 unmeas , and 
a part that depends on our choice of angularities as the observable that arises from soft 
emission into measured jet i, S™ eas (r*). ,5 unmcas i s hence sensitive to the scale A while 
5 mcas (r*) is sensitive to the scale w.;t*. 

In Sec. 6, having calculated all the jet and soft function contributions to 0(a s ), we 
extract the anomalous dimensions and perform renormalization-group (RG) evolution. We 
find the hard anomalous dimension from existing results in the literature. The hard, jet, 
and soft anomalous dimensions have to satisfy the consistency condition Eq. (1.3) in order 
for the physical cross section to be independent of the arbitrary factorization scale \x. Our 
calculations reveal that, as long as the jet separation parameter t Eq. (1.4) between all 
pairs of jets is much larger than 1, the condition is satisfied. 

Even after requiring t > 1, the satisfaction of the consistency condition is non-trivial. 
The hard function knows only about the direction of each jet and the jet function knows 
only the jet size R; the soft function knows about both. Furthermore, it is not sufficient 
only to include regions of phase space where radiation enters the measured jets. We learn 
from our results in this Section that it is crucial to include soft radiation outside of all jets 
with an upper energy cutoff of A. Only after including all of these contributions from the 
various parts of phase space do the jet, hard, and soft anomalous dimensions cancel and 
we arrive at a consistent factorization theorem. 

We conclude Sec. 6 by proposing in Sec. 6.4 a strategy to sum logs due to a hierarchy 
of scales in the soft function, by "refactorizing" it into multiple pieces, each sensitive to a 
single scale, as suggested by the discussion surrounding Fig. 1. Our current implementation 
of this procedure does tame the logarithmic dependence of jet shape distributions on the 
ratio A/uj in our numerical studies, but we leave for further development the resummation 
of all "non-global" logs of ratios of multiple soft scales that begin at NLL and O(a^). 

To help the reader find the results of the calculations in Sec. 4 through Sec. 6 just 
outlined, Table 1 provides a summary with equation numbers. 

In Sec. 7 we compare our resummed perturbative predictions for jet shape distributions 
to the output of a Monte Carlo event generator. We test both the accuracy of these 
predictions and assess the extent to which hadronization corrections affect jet shapes. We 
will illustrate our results in the case of e + e _ — > 3 jets, with the jets constrained to be 
in a configuration where each has equal energy and are maximally separated. In both 
the effective theory and Monte Carlo, we can take the jets to have been produced by an 
underlying hard process e + e~ — > qqg. After placing cuts on jets to ensure each parton 
corresponds to a nearby jet, we measure the angularity jet shape of one of the jets. We 
compare our resummed theoretical predictions with the Monte Carlo output for quark and 
gluon jet shapes with various values of a and R. We find that the dependencies on a and R 
of the shapes of the distribution and the peak value of r a agree well between the theory and 
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Category 


Contribution 


Symbol 


Location 


ATT / \ j_ "1 i • 

INLL) contributions 
before resummation: 


measured quark jet function 


tQ ( \ 


T71 / A 1 1 \ 

Eq. (4.11) 


unmeas. quark jet function 




Eq. (4.17) 


measured giuon jet function 




Eq. (4.25) 


unmeas. gluon jet function 


•Jul 


Eq. (4.26) 


summary of divergent 
parts of soft func. (any t) 





Table 2 


total universal 
soft func. (large t) 


(^•unmeas 


Eq. (5.20) 


total measured 
soft func. (large t) 


S mCaS (T l a ) 


Eq. (5.22) 


dnuniciio Lib dimensions. 






XdUlfci O 


NLO contributions 
after resummation: 


measured jet function 




Eq. (6.42a) 


measured soft function 




Eq. (6.42b) 


unmeas. jet function 




Eq. (6.43a) 


universal soft function 


^unmcas^A) 


Eq. (6.43b) 


Total NLL Distribution: 






Eq. (6.40) 



Table 1: Directory of main results: the fixed-order NLO quark and gluon jet functions for jets 
whose shapes r Q are measured or not; the fixed-order NLO contributions to the soft functions from 
parts of phase space where a soft gluon enters a measured jet, S meas (T a ), or not, S™ meas ; their 
anomalous dimensions; the contributions the jet and soft functions make to the finite part of the 
NLL resummed distributions; and the full NLL resummed jet shape distribution itself. 



Monte Carlo, with small but noticeable corrections due to hadronization. We can estimate 
these corrections by comparing output with hadronization turned on or off in Monte Carlo. 

In Sec. 8, we give our conclusions and outlook. We also collect a number of technical 
details and results for 0{a s ) finite pieces of jet and soft functions in the Appendices. 

Our work is, to our knowledge, the first achieving factorization and resummation of a 
jet observable distribution in an exclusive iV-jet final state defined by a non-hemisphere jet 
algorithm. 4 Having demonstrated the consistency of this factorization for any number of 
quark and gluon jets, measured and unmeasured jets, and phase space cuts in cone and kr- 
type algorithms, and having constructed a framework to resum logarithms of jet shapes in 
the presence of these phase space cuts, we hope to have provided a starting point for future 
precision calculations of many jet observables both in e + e~ and hadron-hadron collisions. 
The case of pp collisions will require a number of modifications, including turning two of 
our outgoing jet functions into incoming "beam functions" introduced in [62]. We leave 
this generalization for future work. 

The reader wishing to follow the general structure of our ideas and logic and understand 
the basis of the final results of the paper without working through all the technical details 
may read Sees. 1 and 2, and then skip to Sec. 7. Some short less technical discussion also 
appears in Sec. 3.4. 

4 Dijet cross sections for cone jets were factorized and resummed in [61]. 
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2. Jet Shapes and Jet Algorithms 

2.1 Jet Shapes 

Event shapes, such as thrust, characterize events based on the distribution of energy in 
the final state by assigning differing weights to events with differing energy distributions. 
Events that are two-jet like, with two very collimated back-to-back jets, produce values of 
the observable at one end of the distribution, while spherical events with a broad energy 
distribution produce values of the observable at the other end of the distribution. While 
event shapes can quantify the global geometry of events, they are not sensitive to the 
detailed structure of jets in the event. Two classes of events may have similar values of 
an event shape but characteristically different structure in terms of number of jets and the 
energy distribution within those jets. 

Jet shapes, which are event shape- like observables applied to single jets, are an effective 
tool to measure the structure of individual jets. These observables can be used to not only 
quantify QCD-like events, but study more complex, non-QCD topologies, as illustrated 
for light quark vs. top quark and Z jets in [12, 60]. Broad jets, with wide-angle energy 
depositions, and very collimated jets, with a narrow energy profile, take on distinct values 
for jet shape observables. In this work, we consider the example of the class of jet shapes 
called angularities, defined in Eq. (1.1) and denoted r . Every value of a corresponds to 
a different jet shape. As a decreases, the angularity weights particles at the periphery 
of the jet more, and is therefore more sensitive to wide-angle radiation. Simultaneous 
measurements of the angularity of a jet for different values of a can be an additional probe 
of the structure of the jet. 

2.2 Jet Algorithms 

A key component of the distribution of jet shapes is the jet algorithm, which builds jets 
from the final state particles in an event. (We are using the term "particle" generically here 
to refer to actual individual tracks, to cells/towers in a calorimeter, to partons in a pertur- 
bative calculation, and to combinations of these objects within a jet.) Since the underlying 
jet is not intrinsically well defined, there is no unique jet algorithm and a wide variety of jet 
algorithms have been proposed and implemented in experiments. The details of each algo- 
rithm are motivated by particular properties desired of jets, and different algorithms have 
different strengths and weaknesses. In this work we will calculate angularity distributions 
for jets coming from a variety of algorithms. Because we calculate (only) at next-to-leading 
order, there are at most 2 particles in a jet, and jet algorithms that implement the same 
phase space cuts at NLO simplify to the same algorithm. At this order the two standard 
classes of algorithms, cone algorithms and recombination algorithms, each simplify to a 
generic jet algorithm at NLO. At NLO the cone algorithms place a constraint on the sep- 
aration between each particle and the jet axis, while the recombination algorithms place a 
constraint on the separation between the two particles. 

Cone algorithms build jets by grouping particles within a fixed geometric shape, the 
cone, and finding "stable" cones. A cone contains all of the particles within an angle R of 
the cone axis, and the angular parameter R sets the size of the jet. In found jets (stable 
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cones), the direction of the total three- momentum of particles in the cone equals the cone 
(jet) axis. Different cone algorithms employ different methods to find stable cones and 
deal with differently the "split /merge" problem of overlapping stable cones. The SISCone 
algorithm [63] is a modern implementation of the cone algorithm that finds all stable cones 
and is free of infrared unsafety issues. In the next-to- leading order calculation we perform, 
there are at most two particles in a jet, and we only consider configurations where all jets 
are well-separated. Therefore, it is straightforward to find all stable cones, there are no 
issues with overlapping stable cones, and the phase space cuts of all cone algorithms are 
equivalent. This simplifies all standard cone algorithms to a generic cone- type algorithm, in 
which each particle is constrained to be within an angle R of the jet axis. For a two-particle 
jet, if we label the particles 1 and 2 and the jet axis n, then the cone-like constraints for 
the two particles to be in a jet are 

cone jet: 9i n < R and #2n < R ■ (2-1) 

This defines our cone-type algorithm. 

Recombination algorithms build jets by recursively merging pairs of particles. Two 
distance metrics, defined by the algorithm, determine when particles are merged and when 
jets are formed. A pairwise metric /O pa ir (the recombination metric) defines a distance 
between pairs of particles, and a single particle metric pj e t (the beam, or promotion, metric) 
defines a distance for each single particle. Using these metrics, a recombination algorithm 
builds jets with the following procedure: 5 

0. Begin with a list L of particles. 

1. Find the smallest distance for all pairs of particles (using p pa i r ) and all single particles 
(using pj ct ). 

2a. If the smallest distance is from a pair, merge those particles together by adding their 
four momenta. Replace the pair in L with the new particle. 

2b. If the smallest distance is from a single particle, promote that particle to a jet and 
remove it from L. 

3. Loop back to step 1 until all particles have been merged into jets. 

The kx, Cambridge- Aachen, and anti-kx algorithms are common recombination algo- 
rithms, and their distance metrics are part of a general class of recombination algorithms. 
For e + e~ colliders, a class of recombination algorithms can be defined by the parameter a: 

p p ^j) = mm(Ef,E?) 6 ^ 

Pie t(i)=E?, (2.2) 

5 This defines an inclusive recombination algorithm more typically applied to hadron-hadron colliders. 
We are applying it here to the simpler case of e + e~ collisions in order to facilitate the eventual transition 
to LHC studies. Exclusive recombination algorithms, more typical of e + e~ collisions, are described along 
with other jet algorithms in [64] and their description in SCET is given in [32]. 
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where a = 1 for kx, for Cambridge- Aachen, and —1 for the anti-kp algorithm. The 
parameter R sets the maximum angle between two particles for a single recombination. 6 
In the multijet configurations we consider the jets are separated by an angle larger than 
R, so that only the pairwise metric is relevant for describing the phase space constraints 
for particles in each jet. For a two-particle (NLO) jet, the only phase space constraint 
imposed by this class of recombination algorithms is that the two particles be separated 
by an angle less than R: 

k T jet: 9 12 <R. (2.3) 

This defines a generic recombination algorithm suitable for our calculation. We will denote 
this as the bp-type algorithm. 

The configurations with two widely separated energetic particles best distinguish cone- 
type jets from kp-type jets at NLO. For instance, the case where the two energetic particles 
are at opposite edges of a cone jet (at an angle 2R apart) is not a single kp jet. However, it 
is important to note that these configurations will not be accurately described in this SCET 
calculation for R 3> A, as such configurations are power suppressed in our description of 
jets. Our concern is in accurately describing the configurations with narrow jets (small r a ), 
and not the wide angle configurations above. 

Because jets are reliable degrees of freedom and provide a useful description of an 
event when they have large energy, in the description of an event we impose a cut A on 
the minimum energy of jets. An A-jet event, therefore, is one where N jets have energy 
larger than the cutoff A, with any number of jets having energy less than the cutoff. In 
our calculation, we impose the same constraint: any jet with energy less than A is not 
considered when we count the number of jets in the final state. This imposes phase space 
cuts: for a gluon radiated outside of all jets in the event, that gluon is required to have 
energy E g < A to maintain the same number of jets in the event. The proper division of 
phase space in calculating the jet and soft functions is a key part of the discussion below, 
and careful treatment of the phase space cuts is needed. 

2.3 Do Jet Algorithms Respect Factorization? 

The factorization theorem places specific requirements on the structure of jet algorithms 
used to describe events. As in Eq. (1.2), the factorization theorem divides the cross section 
into separately calculable hard, jet, and soft functions. The hard function depends only on 
the configuration of jets, while the jet and soft functions describe the degrees of freedom 
in each jet in terms of the observable r. While the soft function is global, the jet function 
depends only on the collinear degrees of freedom in a single jet. The limited dependence 
of the hard and jet functions implies constraints on the jet algorithm. 

Because jets are built from the long distance degrees of freedom arising from evolution 
of energetic partons to lower energies, the configuration of jets in an event depends on 
dynamics across all energy scales. This naively breaks factorization in SCET, since the 

6 We use R for both cone and kT algorithms for ease of notation. For kT, this parameter is sometimes 
referred to as D. We emphasize that having the same size R for different algorithms does not in general 
guarantee the same sized jets. 
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configuration of jets in the hard function would depend on dynamics at low energy in the 
soft function. However, we can describe a jet algorithm that respects factorization, and in 
Sec. 3.4 we will parameterize the power corrections that arise from various algorithms. 

The primary constraint on the jet algorithm in order to satisfy the factorization the- 
orem is that the phase space cuts on the collinear particles in the jet are determined only 
by the collinear degrees of freedom. This essentially ensures that the jet functions are 
independent of dynamics in the soft function. Correspondingly the soft function can only 
know about the jet directions and their color representations. The direction of the jet is 
naturally set by the collinear particles, as soft particles have energy parametrically lower 
than the collinear ones and change the jet direction by a power suppressed amount. The 
further restriction that the phase space cuts on the collinear degrees of freedom are in- 
dependent of the soft degrees of freedom places a strong constraint on the action of the 
jet algorithm. Cone algorithms already implement this constraint: the jet boundary (the 
cone) is determined by the location of the jet axis, which is the direction of the sum of 
all collinear particles up to a power correction. Recombination algorithms, however, are 
constrained by the factorization theorem to operate in a specific way: all collinear particles 
must be recombined before soft particles. As discussed in Sec. 3.4, commonly used algo- 
rithms obey this constraint up to power corrections in the observable for measured jets. 
Of particular note is the anti-kT algorithm, which exhibits behavior very close to what is 
required by the factorization theorem (similarly to the way cone algorithms behave). 

3. Factorization of Jet Shape Distributions in e + e~ to N Jets 

In this Section we formulate a factorization theorem for jet shape distributions in e + e~ 
annihilation to N jets. All the formal aspects we need to describe an A-jet cross section 
appear already in the 3-jet cross section, so we will give explicit details only for that 
case. We will use the framework of Soft-Collinear Effective Theory (SCET), developed in 
[21, 22, 23, 24], to formulate the factorization theorem. We begin with a basic review of 
the relevant aspects of the effective theory. 

3.1 Overview of SCET 

SCET is the effective field theory for QCD with all degrees of freedom integrated out, other 
than those traveling with large energy but small virtuality along a light-like trajectory 
n, and those with small, or soft, momenta in all components. A particularly useful set 
of coordinates is light-cone coordinates, which uses light-like directions n and n, with 
n 2 = h 2 = and n ■ n = 2. In Minkowski coordinates, we take n = (1,0,0,1) and 
n = (1, 0, 0, —1), corresponding to collinear particles moving in the +z direction. A generic 
four- vector p^ can be decomposed into components 

p» = fi- p— + n- p— + p^_. (3.1) 

In terms of these components, p = (n ■ p,n ■ p,p±), collinear and soft momenta scale with 
some small parameter A as 

Pn = E(l,X 2 ,X), p s ~ ^(A 2 , A 2 , A 2 ) , (3.2) 
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where E is a large energy scale, for example, the center-of-mass energy in an e + e~ collision. 
A is then the ratio of the typical transverse momentum of the constituents of the jet to the 
total jet energy. Quark and gluon fields in QCD are divided into collinear and soft effective 
theory fields with these respective momentum scalings: 



l(x) = q n (x) + q s (x) , A*{x) = A»(x) + A£(x) . 



(3.3) 



We factor out a phase containing the largest components of the collinear momentum from 
the fields q n , A n . Defining the "label" momentum pn = n-p n ^ + pj_, where n -p n contains 
the 0(1) part of the large light-cone component of the collinear momentum p n , and p± the 
0(A) transverse component, we can partition the collinear fields q n ,A n into their labeled 
components, 

q n (x) = e-^qnA*) > KV) = E • (3-4) 

p^O p^O 

The sums are over a discrete set of 0(1, A) label momenta into which momentum space is 
partitioned. The bin p = is omitted to avoid double-counting the soft mode in Eq. (3.3) 
[65]. The labeled fields q n , p , A n ^ p now have spacetime fluctuations in x which are conjugate 
to "residual" momenta k of the order EX 2 , describing remaining fluctuations within each 
labeled momentum partition [23, 65]. It will be convenient to define label operators = 
n ■ Vn^/2 +7 3 jf which pick out just the label components of momentum of a collinear field: 



(3.5) 



Ordinary derivatives acting on effective theory fields 4> n ,p(x) are of order EX 2 . 

The final step to construct the effective theory fields is to isolate the two large compo- 
nents of the Dirac spinor q n ^ p for a fermion with lightlike momentum along n. The large 
components £ njP and the small H njP can be separated by the projections 



£n,p — . Qn.p ? ^n,p 



m,p i 



(3.6) 



and we have qn,p — Cn,p ~\~ ^n,p- One can show, substituting these definitions into the QCD 
Lagrangian, that the fields H n>p have an effective mass of order E and can be integrated 
out of the theory. The effective theory Lagrangian at leading order in A is [22, 23, 24] 



^SCET = £■£ + C-A n + &s t 

where the collinear quark Lagrangian is 



in 



in ■ P 



where W n is the Wilson line of collinear gluons, 



W n (x) = ^ eX P 



perms 



n ■ V 



n ■ A n (x) 



(3.7) 



(3.8) 



(3.9) 
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the collinear gluon Lagrangian Ca u is 
1 



Tr 



2g 2 

+ 2Tv{ <■„ 



iV 



+ -TV 

a 



(3.10) 



where c n is the collinear ghost field and a the gauge- fixing parameter; and the soft La- 
grangian C s is 

C s = q s ifl s q s (x) - I Tr G^G 8 Ux) , (3.11) 



which is identical to the form of the full QCD Lagrangian (the usual gauge-fixing terms 
are implicit). In the collinear Lagrangians, we have defined several covariant derivative 
operators, 



D» = d"-igA*-igAi; 



T^ + gA%, W = V^ + in-D 



(3.12) 



In addition, there is an implicit sum over the label momenta of each collinear field and the 
requirement that the total label momentum of each term in the Lagrangian be zero. 

Note the soft quarks do not couple to collinear particles at leading order in A. Mean- 
while, the coupling of the soft gluon field to a collinear field is in the component n-A s only, 
according to Eqs. (3.8) and (3.10), which makes possible the decoupling of such interac- 
tions through a field redefinition of the soft gluon field given in [24]. We will utilize this 
soft-collinear decoupling to simplify the proof of factorization below. 

The SCET Lagrangian Eq. (3.7) may be extended to include collinear particles in 
more than one direction [25]. One adds multiple copies of the collinear quark and gluon 
Lagrangians Eqs. (3.8) and (3.10) together. The collinear fields in each direction rii con- 
stitute their own independent set of quark and gluon fields, and are governed in principle 
by different expansion parameters A associated with the transverse momentum of each jet, 
set either by the angular cut R in the jet algorithm or by the measured value of the jet 
shape T a . Each collinear sector may be paired with its own associated soft field A s with 
momentum of order EX 2 with the appropriate A. For the purposes of keeping the notation 
tractable while proving the factorization theorem in this section, we will for simplicity take 
all A's to be the same, with a single soft gluon field A s coupling to collinear modes in 
all sectors. In Sec. 6.4, we will discuss how to "refactorize" the soft function further into 
separate soft functions each depending only on one of the various possible soft scales. 

The effective theory containing N collinear sectors and the soft sector is appropriate to 
describe QCD processes with strongly-interacting particles collimated in N well-separated 
directions. Thus, in addition to the power counting in the small parameter A within each 
sector, guaranteeing that the particles in each direction are well collimated, we will find in 
calculating an iV-jet cross section the need for another parameter that guarantees that the 
different directions Hi are well separated. This latter condition requires i > 1, where t is 
defined in Eq. (1.4). 7 



7 This condition is a consequence of our insistence on using operators with exactly N directions to create 
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3.2 Jet Shape Distribution in e + e — > 3 Jets 

Consider a 3-jet cross section differential in the jet 3-momenta Pi 2 3> where we measure 
the shape t„ of one of the jets, which we will call jet 1. The full theory cross section for 
e + e _ — >• 7* — >• 3 jets at center-of-mass energy Q is 

dT iX 1A3 = 1 {xl f{0) |0) L/ (2vr) V(Q " px)6n ^ x ^ 

3 (3-13) 
x 5{rl - r a (jet 1)) J] 6 3 {P 3 - P(jet j)) , 

where the J(X) is the jet algorithm acting on final state X, and N{J{X)) is the number 
of jets identified by the algorithm [30]. P(jet j) is the 3-momentum of jet j, and is also a 
function of the output of the jet algorithm J{X). is the leptonic part of the amplitude 
for e + e~ — > 7* — > qqg. The current is 

summing over colors a and flavors /. 

When the three jet directions are well separated, we can match the QCD current j fM (x) 
onto a basis of three-jet operators in SCET [34, 66] . We build these operators from quark 
jet fields Xn, related to collinear quark fields £ n by Xn = Wn£ n , where W n is given by 
Eq. (3.9), and a gluon jet field related to gluons A n by 

B^ = -Wl(V^+Ai)W n . (3.15) 

The matching relation is 

f( x )= E E e ^ 1_p ~ 2+p ~ 3) '*^>i,i^ 

nxn^nz pip2P3 

(3.16) 

with sums over Dirac spinor indices a, f3 and Lorentz index v, and over label directions ^2,3 
and label momenta Pi,2,3- Sums over colors and flavors are implied. We have chosen to 
produce a quark in direction m, antiquark in 712, and gluon in 713. The matching coefficients 
C^p u are found by equating QCD matrix elements of to SCET matrix elements of the 
right-hand side of Eq. (3.16). These coefficients have been found at tree level in [66]. The 
number of independent Dirac and Lorentz structures that can actually appear with nonzero 
coefficients is considerably smaller than suggested by Eq. (3.16) due to symmetries. We 
will keep the form of these coefficients general to make extension to N jets transparent, 
which would require the introduction of a basis of N jet fields in Eq. (3.16), with specified 



the final state. We could move away from the large-t limit and account for corrections to it by using a basis 
of operators with arbitrary numbers of jets and properly accounting for the regions of overlap between an 
iV jet operator and (N ± l)-jet operators. This is outside the scope of the present work, where we limit 
ourselves to kinematics well described by an iV-jet operator, and thus, limit ourselves to the large-t limit. 
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numbers of quark, antiquark, and gluon fields. We will not write the details for an iV-jet 
cross section here, but the procedures are obvious extensions of all the steps in factorizing 
the 3-jet cross section below. 

As a final step before factorization, we redefine the collinear fields to decouple collinear- 
soft interactions in the Lagrangian [24]: 

Xn(x)=Y^(x) X ^(x) (3.17a) 
Xn{x) = ^\x)Y n {x) (3.17b) 
A n (x) = y n (x)A^(x), (3.17c) 

where Y n is a Wilson line of soft gluons along the light-cone direction n, 



Y n (x) = Pexp 



/•oo 

/ dsn ■ A s (ns + x) 
Jo 



(3.18) 



with A s in the fundamental representation. 8 y n is similar but in the adjoint representation. 
The new fields \n\A^ do not have interactions with soft fields in the SCET Lagrangian 
at leading order in A. Henceforth, we use only these redefined fields, but for simplicity 
drop the (0) superscripts. 

The cross section in SCET can now be written, 

dT ij£ 1|2|3 = ^EWW)-^^ - *.Get 1)) ft * 3 (P, - PQet i)) 



x 



X 



E E / dtxe^-n+fr-^C^nw, £1,2,3)^(7^,3; p lj2 , 3 ) 

"1,2,3 Pi, 2, 3 

(0| f{^Ygy^(gB^)T^^(x)} \X) 

{x\ T{ttnK?yZ D (gB££)ffi^&nQ)} 10) • (3-19) 



To proceed to factorize this cross section, it is convenient to rewrite the remaining delta 
functions that depend on the final state X in terms of operators acting on X. Those 
quantities depending on the jet algorithm J can be rewritten in terms of an operator 
containing the momentum flow operator, 

£ M (n)= lim / dtmT^Rn), (3.20) 
Jo 

where T„ v is the energy-momentum tensor, evaluated at time t and position Rn. The 
operator £^(n) measures the flow of four-momentum in the direction n (cf. [29, 68, 69]), 
and the jet algorithm J can be written as an operator J acting on the momentum flow in 
all directions [30]. Correspondingly we can define an operator for the 3- momentum of the 



s The path choice (0 to 00) in Eq. (3.18) is convenient for outgoing particles. The physical cross section 
is independent of whether the path goes to ±00 if the transformation of the external states X is also taken 
into account [67]. 
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jet, P(Jj(i7)). In addition, the event shape r Q (jet 1) can also be expressed as an operator 
T a (Ji(<j)), built from the momentum flow operator, acting on the state \X) (cf. [29]): 

r a (MJ)) = I dve-^-^SrinMi] - ^UW))) ■ (3.21) 



The operator is constructed to count only particles actually entering the jet in direction n\ 
determined by the action of the jet algorithm (for simplicity we will suppress the argument 
J\(j) of f a in the following, but add a superscript for the jet number). Using these 
operators, we can eliminate the X dependence in the delta functions in Eq. (3.19) and 
perform the sum over states X, obtaining 

-2 



° ' ' "1,2,3 pi,2, 3 

3 



i=i 



x T{xT upi Yg yg D (gB%°)Tf g Y^ h xJg^O)} |0) . (3.22) 

The matrix element can be calculated as the sum over cuts of time-ordered Feynman 
graphs, with the delta function operators inside the matrix element enforcing phase space 
restrictions from the jet algorithm and jet shape measurement on the final state created 
by the cut. 

The operators f a and J depend linearly on the energy-momentum tensor, which itself 
splits linearly in SCET into separate collinear and soft pieces, 

t a- = E t ^+ t ^> ( 3 - 23 ) 

i 

which will aid us to factorize the full matrix element in Eq. (3.22) into separate collinear 
and soft matrix elements. To achieve this factorization, however, we must make some more 
approximations : 

1. The contribution of soft particles and residual collinear momenta to the momentum 
P(jet j) of each jet can be neglected, and the jet momentum is just given by the 
label momentum pj of the collinear state \Xj). Thus the energy and jet axis of each 
jet is approximated to be that of the parent collinear parton initiating the jet. In 
particular, the squared mass of the jet is order A 2 compared to its energy. So in 
this approximation we take the jet energy to be equal to the magnitude of its 3- 
momentum. On the other hand, we keep the leading non-zero contribution to the 
angularity even though it is also of order A 2 . These approximations also require that 
we treat the energy of any particles outside all of the jets, and thus the cutoff A, as 
a soft or residual energy. 
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2. The Kronecker delta restricting the total number of jets to 3 can be factored into three 
separate Kronecker deltas restricting the number of jets in each collinear direction m 
to 1, and one Kronecker delta restricting the soft particles not to create an additional 
jet. This approximation requires the separation between jets to be much larger than 
the size of any individual jet so that different jets do not overlap. Factoring the 
restriction on the number of jets in this way is one reason that the parameter tij in 
Eq. (1.4) is required to be large. 

We describe to what extent the algorithms we consider actually satisfy these two approxi- 
mations in Sec. 3.4. For now we assume these approximations and facts hold, which allows 
us to factor the cross section Eq. (3.19), 

drp^23^i23 = W 2 ^ ^ C^K,2,3;o;i l2j 3)C; 7<5p (n 1>2) 3;a; l!2 ,3) 

" ' ' ' ' "1,2,3 ^1,2,3 

X J d 4 xe i(Q-c 1 n 1 /2+ W2 n 2 /2- W 3n3/2). a; f frjfrggtf- _ Tj _ Tg ) 

x <0| X f nU(x)5 N{ j hl d(rj - f^CC) 1°) 5 fa ~ y) ^ ~ n i) (3.24) 

x <0| Xt,-^) S N(J),1XZ-^) 1°) 5 fa - y) s2 ^2 - n 2 ) 

x (0| {gB%X){x)8 N0)tl {gB%* 3 {Q) |0> 6 fa - f ) 6 2 (Sl 3 - n 3 ) 

x (0| Y^Y^T^Y^{x)8 N0)fi 8{r s - r^Y^Y^Y^O) |0> 

We have rewritten the cross section to be differential in E{ (the magnitude of Pj) and Oj 
(the direction of Pj). In the sum over label directions, can be chosen to align with P, 
such that pf = 0. In Eq. (3.24) we have written the label momentum as lo% = rtj • pi. In 
Eq. (1.1) we approximate the jet axis by this rij and the jet energy by n., • pi/2, so that 
they do not depend on soft momenta at all. The operators f™ 1 and are defined to count 
only particles inside the measured jet identified by the algorithm. 

In the soft matrix element in Eq. (3.24), we have introduced the soft Wilson line Y n in 
the antifundamental representation to remove the time- and anti-time-ordering operators 
T,T in Eq. (3.19) [27], and related Wilson lines y n in the adjoint representation to those 
in the fundamental representation by [24] 

yAB T B = Y V T A Yn _ (3 25) 

Defining the jet functions by the relations 



(2k) 



'J ni Mrj,ni ■ fci) (tf) ^ = (0| xRwWn&.AtJ ~ ^X-ZvM l°> 

\ /yet 

(3.26a) 



vX Jn 2 ^(n 2 .k 2 ) [j) S S = (D| X%^ 5 N(J)lXZ*M 1°) ( 3 - 26b ) 
^e^J m (n 3 -k 3 ) 9 ^S AB = -co 3 (0| (gB^S^^gB^g) |0> , (3.26c) 
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and the soft function by 

x Y ni y n3 T B y n3 F n2 (o)|o> 

we can express the cross section Eq. (3.24) as 
da L 2 N F N C C F 



(3.27) 



^1^,2,3^1,2,3 " 6Q 2 

' ' ' ' ".1,2,3 1^1,2,3 



j e i[Q-(wini-wana+W3»i3)/2]-a! (3.28) 

11,2,3 Wl,2,3 

X ^^("1,2,31 Wi,2,3)C* 7 ,5 p (ni,2,3,Wi 1 2,3) fy^ fy) ffl P 
. 3 

x y dT-jdrs^T-i - rj - r s ) n 5 (Ei - y) <5 2 (^ - m) 



i=l 



(2vr) 4 y (2vr) 4 7 (2^) 4 7 (2tt) 

x Jm.WiCrjjrex • fci) J n2 ,-w 2 (n 2 • & 2 ) Jn 3 ,aj 3 («3 • h)S(T S ,r) , 

where now Pj = £j(l,rij). The quark and antiquark jet functions are now for a single 
flavor, and we have summed over Np flavors to obtain the factor in front. The jet functions 
depend only on the smallest component of momentum m ■ ki in each collinear direction. 
The residual and soft momenta appearing in the exponentials can be reabsorbed into the 
label momenta by a series of reparameterizations of the label momenta and directions, 
under which the SCET Lagrangian is invariant [70]. The three-jet operator Eq. (3.16) 
will receive corrections of order A 2 (which we can drop) under the reparameterizations we 
perform below, or can be constructed from the outset to be explicitly reparameterization 
invariant (RPI) [66]. 

First, collect the residual and soft momenta together by defining K = k\ + fc 2 + k% + r. 
We can decompose x in n\ light-cone coordinates, so 

e ~iK-x _ & —i{n\-Kn\-x/2+n\-Kn\-x/2+Kxi-xxi) 29) 

Performing a type- A transformation (in the language of [70]) on the label momentum 

pi = uim/2, 

ui — > oj\ + fi\ • K , (3.30) 
and a type-IB transformation on the vector n\ itself, 

n^m + A 1 , A ± = -—K ± , (3.31) 

shifts the label momentum on the jet function 1 by uj\n\/2 — > {uj\ + fi\ ■ K)n\/2 + K±_\. 
The summation variables n\,oj\ can then be shifted to eliminate n\ ■ K and K±\ from the 
exponentials entirely. We drop all corrections suppressed by A 2 due to these shifts. 

It remains to absorb the n\ ■ K component of residual and soft momentum, appearing 
in the exponential factor e — * n i - ^* i i-x/2_ This cannot be achieved by RPI transformations 
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in the n\ sector since this momentum is purely residual — there is no label momentum 
in this direction. However, in a multijet cross section, h\ can be written as a linear 
combination of n 2 ,ns, and, say, n± 2 (a unit vector transverse to n 2 ,h 2 ), so that RPI 
transformations on u 2 , 0J3 and 77-2 similar to those above can absorb n\ ■ K into the label 
momenta also. Then, the x-dependent residual and soft exponentials all disappear, and we 
can combine the nine superfluous hi ■ ki and k±i integrals with the nine discrete sums over 
label directions and momenta to give continuous integrals over total momenta. Performing 
these with the remaining energy and direction delta functions in Eq. (3.28) and the x 
integral with the remaining exponential gives the momentum conservation delta function 
5 A {Q - E x m - E 2 n 2 - E 3 n 3 ). 

The resulting cross section Eq. (3.28) takes the form 

A 1AV ^ 7 A2 = AT? ^ 1 A2 H ("1.2,35 #1,2,3) / dTjdT S 5(T^ - Tj - T S ) 

dT^dE lj2! 3d z ni t 2,3 dE lj2t3 d 2 n lj2)3 J 

r dni-ki r dn 2 -k 2 f dn 3 -k 3 f d 4 r (3.32) 
X J 2vr J 2tt J 2tt J (2^)4 
x J ni ,2Ei{ T J^i ■ ki)J n2)2 E 2 {n 2 ■ k 2 )J n3)2 E :j (n 3 ■ k 3 )S(r s ,r) , 

where we used that the matching coefficients C^ u (ni;pi) are such that, by construction, 
the right-hand side at tree-level is simply the Born cross section (denoted by a^) for 
e + e~ — > qqg times S(t^). The hard function H = 1 + 0{a s ) is determined by calculating 
the matching coefficients C order-by-order in perturbation theory. 

The above may be easily modified to describe the antiquark or gluon jet angularities, 
by moving the appropriate delta function <5(t* — r a (jet i)) into the J 2 or J 3 jet functions. 
In addition, we may choose from among various jet algorithms. The choice determines 
what G-function restrictions must be inserted into the final state phase space integrations 
created by cutting the jet and soft diagrams to determine which particles make it into the 
jet. As long as the algorithm is such that the approximations enumerated above hold, it 
will not violate factorization of the jet shape cross section. We will discuss factorization 
and its potential breakdown in the context of particular jet algorithms in more detail in 
Sec. 3.4. 

Another check of the validity of the factorization theorem is that the factorized jet and 
soft functions be separately IR safe, which is a stronger condition than the full cross section 
being IR safe. If the observable [54, 71] or algorithm [32] too strongly weights final states 
with narrow jets whose invariant masses are the same as the virtuality of soft particles, 
then the jet and soft functions for the observable in standard SCET with dimensional 
regularization become IR divergent. When this occurs it does not necessarily mean that 
factorization is not possible; but at least not in the standard form derived from the version 
of SCET we utilized above. It could be, for example, that a scheme to further separate 
modes by defining the theory with an explicit cutoff [32] or by factorizing modes by rapidity 
instead of virtuality [65] can restore a version of the factorization theorem. We leave an 
explicit study of which algorithms and observables give IR safe jet and soft functions in 
SCET in dimensional regularization for a separate publication. However, we note here that 
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the algorithms and observables (r a for a < 1) that we consider in this paper, at least at 
NLO, do give rise to IR safe jet and soft functions. 



3.3 Jet Shapes in e + e — > N jets 

To generalize the result Eq. (3.32) to an arbitary number N of jets, we simply add more 
quark and gluon jet fields to the operator matching in Eq. (3.16), resulting in the corre- 
sponding number of additional quark and gluon jet functions in Eq. (3.32), along with a 
hard coefficient and a soft function for N jet directions. Consider an event with 2N q quark 
and antiquark jets and N g gluon jets, where 2N q + N g = N . Furthermore, we can choose 
to measure the angularity shape for any number of these jets. Achieving a factorization 
theorem that remains consistent for any of these combinations is a nontrivial task and thus 
a powerful test of the validity of the effective theory. 

For an N = 2N q + N g jet event, we generalize the matching of the QCD current in 
Eq. (3.16) to: 



? -M(V) - V V e i(pi+-+PN)-x C ^- a N q bi-bN q A 1 -A Ng 



ni— njv Pi—pn 



(3.33) 



i=l j=l k=l 



with implicit sums over Dirac spinor indices cti, /3j, Lorentz indices u^, (anti-)fundamental 
color indices a^, bj, and adjoint color indices A^. The N jet cross section differential in M 
jet shapes, with M < N, factors in SCET into the form 



dajEuni] ■ ■ ■ ;E N ,n N ) 
drl ■ ■ ■ dr™ 

ai aM 



<7 (0, (£i,ni; ■■■;E N , n N )H a ^(n 1 ,E 1 ; ■■■n N ,E N ) 

M 



X 

1=1 

x 4\,2 E yj^i ■*!)■■• 4%,2EjTf;n M ■ k M ) 
jfu+i 

d 4 r 



driN ■ k]\[ 
2^ 



x / ^^^(^....^r), (3.34) 

where <t(°) is the Born cross section for e + e to N q quarks, N q antiquarks, and N g gluons; 
the color indices on the hard and soft functions H and S allow for color mixing; and fi is the 
flavor of each jet function (quark, antiquark, or gluon). H itself is determined by calculating 
the matching coefficients C in Eq. (3.33). The jet functions have the same definitions given 
in Eq. (3.26), and the soft function is given by the appropriate generalization of Eq. (3.27) 
with Wilson lines in the directions and color representations corresponding to the choice 
of fields in Eq. (3.33). We rearrange the order of flavor and color indices in the hard and 
soft functions to agree with the choices of flavor indices on the jet functions. 
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3.4 Do Jet Algorithms Induce Large Power Corrections to Factorization? 

In this section we explore when power corrections to the factorization theorem above be- 
come large, in particular those that are due to the action of jet algorithms. We will argue 
that power corrections to jet angularities induced by the commonly-used cone and recom- 
bination algorithms remain suppressed as long as R is sufficiently large. In particular, we 
need in general that R satisfies R > A and, for the case of the kx algorithm, we require 
R > A. 9 These power corrections are associated with assumptions we made about the 
action of the jet algorithm on final states in deriving Eq. (3.34). In general, the size of 
these power corrections depends both on the algorithm and the observable. 

Power corrections to the px of a jet arising from perturbative emissions (as well as 
from hadronization and the underlying event in pp collisions) for various jet algorithms 
were explored in [36]. These power corrections arise for similar reasons as those we discuss 
below, namely, perturbative emissions changing which partons get combined into the jet. 
Ref. [36] finds that such power corrections scale like lni? for small R. This result is 
consistent with our qualitative discussion below, where we argue that power corrections to 
jet angularities arising from the jet algorithms we use are minimized when R is sufficiently 
large. For us, R should be at least 0(A) (and for the case of the kx algorithm, R 3> A). 

One set of power corrections that is independent of the choice of algorithm arises 
from the approximation of setting the jet axis equal to the label direction n. Since this 
neglects the effects of soft particles, it is valid up to 0(A 2 ) corrections. It was argued in 
Refs. [17, 29, 41] for the case of hemisphere jet algorithms that these corrections in turn 
induce corrections to the angularity r a of order A 2 ( 2_a ), which, for r a ~ A 2 , are subleading 
for a < 1. Essentially the same arguments can be applied to all of the algorithms we 
consider. 

Jet algorithm-dependent power corrections arise from the difference in soft particles 
included in a jet by a given algorithm and those included in the soft function in the 
factorization theorem. The algorithms also differ amongst themselves in which soft particles 
they include in a jet. For observables that scale as 0(1), such as the jet energies and 3- 
momenta, the contribution of soft momenta can be neglected since they scale as 0(A 2 ). 
Clearly then, these observables are not dependent on our choice of jet algorithm and so 
the assumptions we made about factorization of the algorithm in deriving Eq. (3.34) are 
trivially satisfied. 

However, for observables that scale as 0(A 2 ) such as angularities, soft contributions 
become important and so the details of the algorithms we consider become relevant. We 
now estimate how closely the phase space region included in the soft function in the SCET 
factorization theorem approximates the region of soft particles actually included by the 
jet algorithm. We will argue that unless R > 0(A) for the anti-kp and infrared-safe cone 
algorithms, and R 3> A for the kx algorithm, the mismatch in areas is sufficiently large 

9 We noted earlier that there may be different A's for the SCET modes describing different jets. For 
measured jets, A 2 ~ r , while for unmeasured jets, A ~ tan(i?/2), and for soft gluons outside jets, A 2 ~ A/Q. 
In this section, we mean by A the expansion parameter associated with measured jets, and ensuring R is 
much larger than this A. But if R is too large, the separation parameter t oc 1/ tan(i?/2) becomes too small. 
We will consider R ~ 0.4 to 1 to be safe. 
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[^| Soft Region Common to Both 

Algorithm and SCET Factorization Thai. 

Additional Soft Region of Algorithm 
X Parent Location (before collinear splitting) 
# Daughter Location (after collinear splitting) 



Figure 2: Difference between regions of soft radiation included in the SCET factorization theorem 
and the actual (A) ItT and (B) anti-kr algorithms. We illustrate how the regions of soft radiation 
included by the algorithms change when a single, energetic parent particle splits into two collinear 
daughters. Both the algorithm and the soft function merge soft particles contained in the large 
white circle. The algorithm also merges the hatched area and hence contains a region of phase 
space which is different than that included in the SCET soft function (since, as explained in the 
text, the shape and size of the region used in SCET cannot depend on the details of collinear 
splittings). 

so as to cause a leading-order power correction to the measured jet angularities. For R 
larger than these bounds, the corrections are negligible. This miscounting arises due to the 
fact that factorization requires that collinear particles be combined first, and that the soft 
function only knows about the parent collinear direction. When algorithms do not obey 
this ordering, factorization may be violated. 

To determine the size of the soft particle phase space region for each jet algorithm that 
is not included by the factorization theorem, we consider the situation depicted in Fig. 2. 
A parent collinear particle splits into two daughter collinear particles. In the factorization 
theorem, since collinear particles are combined first, the region of phase space where soft 
particles are combined into the jet is a circle of radius R about the parent particle direction. 
However, in a jet algorithm, soft particles in an additional region outside of this may also 
be combined into the jet (the hatched regions in Fig. 2). If the area of this region is of the 
same order as the area included by the factorization theorem, then the power corrections 
to jet angularities induced by the jet algorithm will be leading-order. 

Because soft particles have momenta that are parametrically smaller than collinear 
particle momenta, we determine the omitted region of soft particle phase space by con- 
sidering the dominant action of the jet algorithm. The kr algorithm serves as a useful 
example. The kx metric between a pair of soft particles is O(X 2 )0 SS , the metric between a 
soft particle and collinear particle is 0(X 2 )9 CS , and the metric between a pair of collinear 
particles is O(X )6 cc . Therefore, collinear-collinear recombinations only occur if the angle 
9 CC between the collinear particles is smaller than the separation between any soft parti- 
cle and its nearest neighbor by a factor of 0(X 2 ). Given that the typical angle between 
collinear particles is 0(A), the dominant action of the jet algorithm is to first merge all 
soft particles with their nearest neighbors, while collinear-collinear recombinations occur 
late in the operation of the algorithm. This description will suffice to accurately determine 
the area of the omitted soft phase space for the kx algorithm. Since collinear particles 
are combined last, on average soft particles within circles of radius R about the collinear 
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daughter particles are included in the kx algorithm jet, as shown in Fig. 2A. The area 
that the kx algorithm includes that the soft function does not is represented by the hatched 
region, which is an area of O(0ijR). This area must be parametrically smaller than the 
area included by the soft function (of 0(R 2 )) for the associated power corrections to be 
small. We thus require that R 3> Oij ~ A in the SCET power counting. 

The anti-kx algorithm combines particles in a manner that is closer to respecting 
factorization. It finds the hardest particle first and merges particles at successively larger 
distances from this particle. For the example of two collinear daughters, it will merge all 
soft particles with the hardest daughter that are closer than the distance to the softer 
daughter before merging the two daughters and then merging all soft particles a distance 
R from the merged daughters (i.e., the parent particle), as shown in Fig. 2B. As the Figure 
illustrates, the hatched area of the anti-kx jet tends to be smaller than that of the kx 
jet. In fact, for R > 30y/2, this region vanishes completely (and this case of having only 
two collinear daughters is a worst-case scenario). This leads us to expect that, for any 
number of collinear splittings, for R > A (i.e., not necessarily parametrically larger), power 
corrections due the action of the anti-kx algorithm vanish. 

Cone algorithms such as the SIS Cone algorithm can also include regions that differ 
from the lowest-order region at higher orders in perturbation theory. We now argue that 
an arithmetic bound R > A is sufficient to minimize the power correction from these 
differences, as for the anti-kx algorithm. This situation arises due to the fact that stable 
solutions may exist with overlapping cones when collinear splittings are larger than the cone 
radius, i.e., R < 9ij. In these cases, the amount of radiation that falls into the overlapping 
region is used to decide whether the cones are split or merged. In either case, the boundary 
of the resulting jet(s) has roughly the appearance of Fig. 2A and the difference between 
the region of soft radiation assumed in SCET and that by the actual algorithm is 0(1). 
However, for R > 9{j ~ A for the case of a single collinear splitting, all of the collinear 
radiation lies within a region of size R and there will always be a stable cone that includes 
this radiation and thus the algorithm and the SCET soft function will be sensitive to soft 
particles in the same region of phase space. 

In summary, we have argued that for all the algorithms we consider (kx, anti-kx, 
infrared-safe cone), power corrections are negligible for sufficiently large R. While anti-kx 
and cone allow simply R > A instead of R 3> A as for kx, we will in fact always consider 
R S> A (for the A in a measured jet sector) in the remainder of this paper, guaranteeing 
small power corrections for all these algorithms. (R still determines the scale A in an 
unmeasured jet sector.) Our focus will remain on resumming logs of jet shapes such as 
angularities r a in the presence of jet algorithms, without worrying about resumming logs 
of R themselves. 11 



1 Soft particles in this region can also be removed from this region by merging with other soft particles 
outside of the region and vice- versa, but this average area suffices for our discussion. 

11 Because small R (< 0.3) jets cannot be well resolved in current experiments, resummation of logarithms 
of R is not of primary practical importance in the near future. 
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4. Jet Functions at 0(a s ) for Jet Shapes 

In this section, we calculate the quark and gluon jet functions for jet shapes at next-to- 
leading order in perturbation theory. The jet functions can be divided into two categories: 
those for measured jets, which are fixed to have a specific angularity r a , and those for 
unmeasured jets, which are not. We will denote the quark jet function by J<J, the gluon 
jet function by J2, where u is the label momentum, and a jet function J g ' 9 (r a ) with an 
argument of r a denotes a measured jet. We will calculate the jet functions for the two 
classes of jet algorithms, kp-type and cone-type algorithms. 

In the course of these calculations, we will demonstrate the crucial role of zero-bin 
subtractions [65] from collinear jet functions in obtaining the consistent anomalous dimen- 
sions and the correct finite parts. In this case zero-bin subtractions are not merely scaleless 
integrals converting IR to UV divergences, but in fact contribute part (sometimes the most 
important part) of the correct nonzero result, as was already pointed out by [32, 72]. The 
relation of zero-bin subtractions in SCET to eikonal jet subtractions from soft functions in 
traditional methods of QCD factorization was explored in [41, 73, 74]. In addition, we find 
that the zero-bin subtraction removes the contribution of collinear emissions that escape a 
jet, leaving only power-suppressed pieces in A/wj. 

4.1 Phase Space Cuts 

To calculate the jet functions for a particular algorithm, 
we must impose phase space restrictions in the matrix ele- 
ment. From the jet function definitions, Eq. (3.26), these cuts 
take two forms. One kind, imposed by the operator $ N (j^ 1 
in Eq. (3.26), is common to every jet function. It is the set 
of phase space restrictions related to the jet algorithm, and 
requires exactly one jet to arise from each collinear sector of 
SCET. The other, imposed by the operator S(r a — r a ), is im- 
plemented only on measured jets and restricts the kinematics 
of the cut final states to produce a fixed value of the jet shape. 
In this section we describe these phase space cuts in detail. 

The typical form of the NLO diagrams in the jet functions 
is shown in Fig. 3. As shown in the figure, the momentum 
flowing through the graph has label momentum l~ = n ■ I = oj and residual momentum 
Z + = n ■ I, and the loop momentum is q. We will label "particle 1" as the particle in the 
loop with momentum q and "particle 2" as the particle in the loop with momentum I — q. 
For the quark jet, we take particle 1 as the emitted gluon and particle 2 as the quark. 

As usual, the total forward scattering matrix element can be written as a sum over 
all cuts. Cutting through the loops corresponds to the interference of two real emission 
diagrams, each with two final state particles, whereas cutting through a lone propagator 
that is connected to a current corresponds to the interference between a tree-level diagram 
and a virtual diagram, each with a single final state particle. Thus, the phase space 
restrictions and measurements we impose act differently depending on where the diagrams 




Figure 3: A representa- 
tive diagram for the NLO 
quark and gluon jet func- 
tions. The incoming mo- 
mentum is I = T^uj+Tjl + and 
particles in the loop carry 
momentum q ( "particle 1" ) 
and I — q ( "particle 2" ) . 
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are cut. In addition, since we will be working in dimensional regularization (with d = 4— 2e), 
which sets scaleless integrals to zero, the only diagrams that contribute are the cuts through 
the loops. This means that we only need to focus on the form of phase-space restrictions 
and angularities in the case of final states with two particles. 

The regions of phase space for two particles created by cutting through a loop in the 
jet function diagrams can be divided into three contributions: 

1. Both particles are inside the jet. 

2. Particle 1 exits the jet with energy E\ < A. 

3. Particle 2 exits the jet with energy Ei < A. 

In contributions (2) and (3), the jet has only one particle, which is the remaining particle 
with E > A. 

It is well known that collinear integrations of jet functions can be allowed to extend 
over all values of loop momenta so long as a "zero-bin subtraction" is taken from the result 
to avoid double counting the soft region already accounted for in the soft function [65] . We 
will demonstrate that contributions (2) and (3) are power suppressed by 0(A/ui), which 
scales as A 2 , after the zero-bin subtraction. 

The phase space cuts that enforce both particles to be in the jet depend on the jet 
algorithm. There are two classes of jet algorithms that we consider, cone-type algorithms 
and (inclusive) kr-type algorithms, and all the algorithms in each class yield the same 
phase space cuts. We label the phase space restrictions as CO ne and 0k x , generically a i g - 
For the cone- type algorithms, 

e conc = 9 cone (g, /+) = 9 (tan 2 f > Jp) @ (*an 2 | > • (4-1) 

These functions demand that both particles are within R of the label direction. For the 
kp-type algorithms, the only restriction is that the relative angle of the particles be less 
than R: 

a- 1- a 2 



©k T = ®k T (qJ + ) = ( cosR < 

q^l 2 + q 2 -2q-l / 

= e(wf > ). (4.2) 

In the second line we took the collinear scaling of q (q + Cg"). While this is not strictly 
needed, it makes the calculations significantly simpler. 

For the phase space restrictions of zero-bin subtractions, we take the soft limit of 
the above restrictions. The zero-bin subtractions are the same for all the algorithms we 
consider. For the case of particle 1, which has momentum q, the zero-bin phase space cuts 
are given by 

ei?^ = 0i° o L = e£ = 0(tan 2 f (4.3) 
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Figure 4: Diagrams contributing to the quark jet function. (A) and (B) Wilson line emission 
diagrams; (C) and (D) QCD-like diagrams. The momentum assignments are the same as in Fig. 3. 

The zero bin of particle 2 is given by the replacement q — >■ I — q. 

For all the jet algorithms we consider, the zero-bin subtractions of the unmeasured jet 
functions are scaleless integrals. 12 However, for the measured jet functions, the zero-bin 
subtractions give nonzero contributions that are needed for the consistency of the effective 
theory. 

In the case of a measured jet, in addition to the phase space restrictions we also demand 
that the jet contributes to the angularity by an amount r a with the use of the delta function 
<5f? . = S{r a — f a ), which is given in terms of q and I by 

SR - 5 R (q,l+) =s(r a - i( W - q-T'\l + - q^ 2 - V) a/ V) W2 ) • (4-4) 

In the zero-bin subtraction of particle 1, the on-shell conditions can be used to write the 
corresponding zero-bin <5-function as 

4 0) = * (ra ~ ^-T l % + ) 1 - a/2 ) , (4-5) 
(and for particle 2 with q — >■ I — q). 

4.2 Quark Jet Function 

The diagrams corresponding to the quark jet function are shown in Fig. 4. The fully 
inclusive quark jet function is defined as 

J dSe*' <0| O*)x£,(0) 10) = 5 ab (J )^ J q .(l + ) , (4-6) 

and has been computed to NLO (see, e.g., [75, 76]) and to NNLO [77]. Below we compute 
the quark jet function at NLO with phase space cuts for the jet algorithm for both the 
measured jet, J^(r a ), and the unmeasured jet, J%. As discussed above, we will find that 
the only nonzero contributions come from cuts through the loop when both cut particles 
are inside the jet. 

12 Note that algorithms do exist that give nonzero zero-bin contributions to unmeasured jet functions [32]. 
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4.2.1 Measured Quark Jet 

The measured quark jet function includes contributions from naive Wilson line graphs (A) 
and (B) and QCD-like graphs (C) and (D) in Fig. 4. The sum of these contributions is 



J q ul (r a )=g 2 ^C F 



dl+ 1 



d d q 
(2vr) d 



x Q(q~)Q(q + )2TT5 (l 



l + l + 
q lo 



- Q 

)e(i+ 



2n5(q-q + - q\) 



q + ) ©alg^/J 



(4.7) 



The contribution proportional to d — 2 comes from the QCD-like graphs (C) and (D) in 
Fig. 4. Only the Wilson line graphs have a nonzero zero-bin limit, which comes from taking 
the scaling limit q ~ A 2 of the naive contribution: 

d d q 1 



dr l 



(2?r) d q- 



: 27T5(q-q + -q 2 ± )e(q-)Q(q + ) 



2tt5 (Z+ - ,+) 9(i + - q + ) eM> 



(4.8) 



All jet algorithms that we use yield the same zero-bin contribution, since the phase space 
cuts are the same. 

To evaluate these integrals, we can analytically extract the coefficient of 5(r a ) by 
integrating over r a and using the fact that the remainder is a plus distribution, as defined 
in Eq. (A. 2). We find the naive contribution is 



o s Cf 
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47T/i 2 
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? + 2e 



(4.9) 



The only difference between the jet algorithms that we consider resides in the finite distri- 
bution J^i g (T a ), which is a complicated function of r a that we give in Appendix A. Note 
that the divergent part of the naive contribution is proportional to S(r a ). This is due to 
the fact that the jet algorithm regulates the distribution for r a > 0. The divergent plus 
distributions come entirely from the zero-bin subtraction, which is given by 
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7T 



ra 



(l-a)er a 1+2e 



(4.10) 



Adding the leading-order contribution to all of the NLO graphs and expanding in 
powers of e, adopting the MS scheme, we find the total quark jet function, 



J q M = S(r a ) + J q M ~ Jf\r a ) = 1 + 



OsCf 



OsCf 
TT 



1 



1 1 



1 - a e a ' 1 
1 6(r a ) 



- 1 ii 2 s 



+ 



a e 



4r 



>5{T a ) 



e 1 



a r„ 



^4l s (r a ). (4.11) 



This agrees with the standard jet function J(k + ) given in [75, 76] by setting a = and 
k + = 0JT a . We have shown the divergent terms explicitly, and collect the finite pieces in 
Jgjg(T a ), which we give in Eq. (A. 14). Note that there is no jet algorithm dependence in 
the divergent parts of the jet function at this order in perturbation theory. 
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4.2.2 Gluon Outside Measured Quark Jet 

In this section we calculate the contribution to the quark jet function from the region of 
phase space in which the gluon exits the jet carrying an energy E g < A. This cut causes 
the contribution to be power suppressed by A/ui, which scales as A 2 . However, we elect 
to evaluate this case explicitly as it provides a clear example of the zero-bin subtraction 
giving the proper scaling to the total contribution. We only evaluate this contribution 
for the cone algorithm; the details of the kx algorithm calculation are similar. Note that 
the contribution when the quark is out of the jet is power suppressed at the level of the 
Lagrangian given in Sec. 3.1, in which soft quarks do not couple to collinear partons at 
leading order in A. 

For the cone algorithm, the gluon exits the jet when the angle between the jet axis, 
ni, and the gluon is greater than R. When the gluon is not in the jet, the cone axis is the 
quark direction, and so it makes no contribution to the angularity. Therefore, this region 
of phase space contributes only to the <5(r a ) part of the angularity distribution. 

For the naive contributions, requiring the gluon to be outside the jet and have energy 
less than A, we have the integral 



JSTIt.) = 9 VC F J J ^ ^4- + (4 - 2)—^ j 2^( 9 -,+ - 9l 

x e(q-)e{q + )2ir5 (l + -q+ Q(u - q~)@(l + - q + ) 

V w-g J 

x e(^-tan 2 !^ (2A - g") *(r ) . (4.12) 

Note that the theta function requiring q~ < 2A is more restrictive than q~ < oj. Evaluating 
Eq. (4.12) yields a contribution that scales with A only below the leading term in 1/e: 

rout a s C F 1 ( 4tt// 2 V t/ _^1 1 ^4A 2A 2 ^ 8A^ 



The zero-bin subtraction of Eq. 4.12 is 

zr*>M = / / (£+(<<- » + - A) 

x e(q-)Q(q + )2ir5 (l+ - q+) 9 - tan 2 ^) 9 (2A - q~) 5(r a ) . (4.14) 
Evaluating Eq. (4.14), we find the zero bin will exactly remove the leading term in 1/e: 

Jr ut(0) (r a ) = - as ° F , 1 , ( ^\ \s(r a )\ (4.15) 
K J 2vr r(l-e) l v (2Atanf) 2 y / K e 2 V ' 

Therefore, the difference is power suppressed only after the zero bin is included. Because 
other contributions when one particle is outside of the jet are similarly power suppressed, 
we will drop them in our remaining discussion of the jet functions. 
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4.2.3 Unmeasured Quark Jet 

When the angularity of a jet is not measured, the jet function has no r a dependence. The 
naive and zero-bin contributions are the same as Eqs. (4.7) and (4.8) except for the factor 
of 5r. The zero-bin contribution is 



x 2it5 (l + - q+) Q(l + - q + ) 



(4.16) 



alg • 

This integral is scaleless and therefore equal to in dimensional regularization. This implies 
that the NLO part of the quark jet function for an unmeasured jet is just the naive result. 
We find, making the divergent part explicit, in the MS scheme, 

where the finite part J^ lg is given in Eq. (A. 18). 13 
4.3 Gluon Jet Function 

The diagrams needed for the gluon jet function at NLO are shown in Fig. 5. The fully 
inclusive jet function, defined as 

J #xe*- <0| Bfi,(x)B£,(p) |0> = -±g1 v 6 AB J°(l+) , (4.18) 

(with JSj{1 + ) normalized to 2ir5(l + ) at tree-level) has been calculated to NLO in Feynman 
gauge in [34, 78, 79] and was reported to give the same result in both and light-cone 
gauges in [35] . Since our phase space restrictions and the observables act differently on cuts 
through loops and on cuts through external propagators, it is useful to reproduce these 
results by explicitly cutting the diagrams. 

After some algebra, we find that the sum of all cuts through the loops of the naive 
collinear graphs gives 



/ 



dl+ Tari+\ - „2e 2 9 2 f dl 



1 \ 1 — e ujI + J \ q oj - q uol + J 



(4.19) 



We also record the zero bin that needs to be subtracted from Eq. (4.19). To leading-power 
it is given by 



^^-^(^r/fr/A 



5(q 2 Ml + -q + )e(q-)±: 



+ S((l-q) 2 )5(q + )Q(u- 



(4.20) 



13 The unmeasured jet function Eq. (4.17) is not simply obtained by integrating the measured jet function 
Eq. (4.11) over r a . This is due to the different relative scaling of R with the SCET expansion parameter 
\i in a measured and unmeasured jet sector, as noted earlier. Namely, R ~ in a measured jet sector 
(where A ~ y/ra) while A^ ~ tan(i?/2) in an unmeasured jet sector. 
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Figure 5: Diagrams contributing to the gluon jet function. (A) sunset and (B) tadpole gluon 
loops; (C) ghost loop; (D) sunset and (E) tadpole collinear quark loops; (F) and (G) Wilson line 
emission loops. Diagrams (F) and (G) each have mirror diagrams (not shown). The momentum 
assignments are the same as in Fig. 3. 

Without inserting any additional constraints, this integral is scaleless and zero in dimen- 
sional regularization. Therefore, in the absence of phase-space restrictions, the naive inte- 
gral Eq. (4.19) gives the standard (inclusive) gluon jet function 
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v-l-e 



is + f e 



(4.21) 



in the MS scheme. The measured and unmeasured jet functions are obtained by inserting 
©aig^R and ©aig> respectively, into Eqs. (4.19) and (4.20). 

4.3.1 Measured Gluon Jet 

The naive contribution to the measured gluon jet can be written as 



J 9 M 



a s 1 / 47Tfi 2 \ e 1 / 1 



2vrr(l - e) V J 1 - | \r a 
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@ai g (a:) , 



where x = q /co. This gives 

J9M = 2iw^7) Lw#) 6{Ta) 



+ 2vr ^alg( r a) 



(4.23) 



where, as for the quark jet function, the finite distributions J^ s (T a ) differ among the 
algorithms we consider. They are given in Appendix A. 
The zero-bin result is 



a s C A 1 AVtan^f V/^ 1 
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l+2e 



1 



(4.24) 
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Subtracting the zero-bin from the naive integral and adding the leading-order contribution, 
we obtain in MS 



J r i(r a ) = 5{r a ) + J^r a )-J^){r a ) 
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(4.25) 



7T 1 — a e \ r a / , 27T 
The finite distribution J^ g {j a ) is given in Eq. (A. 14). 
4.3.2 Unmeasured Gluon Jet 

As for the quark jet function, for unmeasured jets the naive and zero-bin contributions are 
given by the measured jet contributions without the 5r constraint. Also, as it was for the 
quark jet function, the zero-bin contribution to the unmeasured jet function is a scaleless 
integral. Thus, the final answer is just the naive result, which is given by 



Jl = 1 + 
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(4.26) 



with the finite part given in Eq. (A. 29) in the Appendix. 



5. Soft Functions at 0(a s ) for Jet Shapes 

In this section we compute the soft function for a generic jet event. In Sec. 5.1, we 
describe the phase space cuts that we impose on soft particles emitted into the final state. 
We then give an outline of how we disentangle the various contributions to the iV-jet 
soft function in Sec. 5.2 and proceed to calculate these contributions in the remaining 
subsections. 



5.1 Phase Space Cuts 

Soft particles in the final state must satisfy the phase space cuts required by the operator 
^N(j) o m (^•^'' 7 )i wri ich constrains the soft particles to not create an extra jet. A soft 
particle is allowed in the final state if it is emitted into one of the jets as defined by the 
jet algorithm, or if it is not in a jet but has energy less than a cutoff A. At NLO, only a 
single soft gluon can be emitted. Therefore, for both cone-type and kp-type algorithms, 
the constraint that the soft gluon be in a jet is simply that the angle of the gluon with 
respect to the jet axis satisfies 9 g j < R. Thus, our requirement on soft gluons is that they 
obey one of the two following conditions: 

in jet i: 9 gJi < R 
out of all jets: E g < A and 9 g j i > R for all % . (5-1) 
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(A) (B) (C) (D) 

Figure 6: Soft function real-emission diagrams. Diagrams (A) and (C) are interference diagrams 
of Wilson line emission from lines i and j and (B) and (D) are from lines i and k. The shaded 
region in the center represents the region of phase space corresponding to jet k defined by the jet 
algorithm, and so the gluons in diagrams (A) and (B) are inside jet k and those in (C) and (D) are 
not. Each diagram has a corresponding mirror diagram (not shown). 

These conditions can be written in terms of theta functions on the gluon momentum k. 
We denote the energy restriction for out-of-jet gluons as 

9 A = Q(k° < A) , (5.2) 

and we denote the requirement that a gluon be in jet i in terms of the light-cone components 
k^ about the direction of jet i, TTj . clS 

&R = 6 (jp < tan 2 |) . (5.3) 

For the case when the soft gluon is in a measured jet, we demand that it contributes 
an amount r a to the angularity of a jet with label momentum uj with the use of the delta 
function 

Sr = s(r a - \k-) a l\k + ) l - a l 2 ^ . (5.4) 
5.2 Calculation of contributions to the A^-Jet Soft Function 

The topology of the various graphs that we need to compute the soft function is shown 
in Fig. 6. The next-to-leading order part Sn\ of the soft function S is the sum of the 
interference of soft gluon emissions from Wilson lines in directions i and j, Sij, over all 
lines i and j with i ^ j (since for i = j, the diagram is proportional to nf = 0), 

It is conceptually straightforward to see that Sij is just the sum of the following three 
classifications of the final state cut gluon: 

• The gluon is in a measured jet and thus contributes to the jet angularity. 

• The gluon is outside of all the jets and has energy E g < A. 

• The gluon is in an unmeasured jet and has any energy. 
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However, the second contribution is technically difficult to compute due to the complicated 
shape of the space with all jets cut out of it, like Swiss cheese. A division of phase space 
leading to a simpler calculation is the following: 

• 5^ eas (r a fc ): The gluon is in measured jet k and contributes to the jet's angularity . 

• Sij'. The gluon is in jet k with energy E g > A (and does not contribute to r„). 

• Sij'. The gluon is in jet k with energy E g < A (and does not contribute to r„). 

• Slj cl : The gluon is anywhere with E g < A (and does not contribute to any angularity). 

In terms of these pieces, the NLO soft function with M measured jets and N — M unmea- 
sured jets is given by 
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(5.6) 



From the definitions above, it is easy to see that the term in large parentheses on the 
second line is equivalent to the sum of the last two contributions on the original list above, 
i.e., the contributions from a gluon not in any jet with E g < A and from a gluon in an 
unmeasured jet with any energy. 

We can simplify this expression by noting that the contribution from a gluon in jet 
k with no energy restriction involves a scaleless integral over the energy that vanishes in 
dimensional regularization and thus 
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Using this relation, the soft function simplifies to 
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(5.8) 



where the first term in Eq. (5.8) is a universal contribution that is needed for every iV-jet 
observable, defined as 
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The second term, defined as, 
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depends on our choice of angularities as the observable. 

We now proceed to set up the one-loop expressions for the contributions in Eq. (5.8). 
The integrals involved are highly nontrivial and so in this section we simply quote the result 
of each integral, referring the reader to Appendix B for details. Most of these integrals are 
most easily written in terms of the variable ty, defined in Eq. (1.4) as tij = tan / tan 5, 
where ipij is the angle between jet directions i and j. (That is, rij • rij = 1 — cosipij.) In 
Table 2, we summarize the divergent parts of the soft function. 

The Feynman rules for the emission of a soft gluon from fundamental- and adjoint- 
representation Wilson lines (corresponding to quark and gluon jets, respectively) have 
the same kinematic structure. The difference is entirely encoded in the color charge of the 
Wilson lines which, using the color space formalism of [80, 81], we denote as Tj for emission 
from Wilson line i. Thus, we can consider the iV-jet soft function without specifying the 
color representation of the final-state jets. 



5.2.1 Inclusive Contribution: S^ cl 

The contribution to the soft function from a gluon going in any direction with energy 
Eg < A is given by the integral 

Sif = - 9 V< T, . T, / ^ (ni ^. fe) ^W) 6 A . (5.11) 
We evaulate this integral in Sec. B.l of the Appendix and find 

* -^m w (? - \ - ^ - t - * - s^)) • <-) 

Note that this calculation is related to the inclusive, timelike soft function that has 
applications elsewhere (see, e.g., [82, 83, 84]), generalized for non back-to-back jets: 

- J„2e T in f d d k rii ■ rij n ^xfi.2\aifL.0\xri,0 



dA 



5.2.2 Soft gluon inside jet k with E g > A: 

Using Eq. (5.7), the contribution from a gluon emitted into jet k from lines i and j is 
given by the integral 

^v-t.-t,/ ^ ( ^ )( r;v ^ 2)e( * 0)9 ^- (5 - i4) 

Much like for the Sfj cas contribution, if k = there is an additional divergence (arising 
from rife • k — > 0) relative to the case k ^ and so we evaluate these two cases separately 
below. 
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Case 1: k = i or j The integrals for this case are performed in Sec. B.2 of the Appendix, 
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Case 2: k ^ i,j These contributions are at most 1/e divergent since the matrix element 
does not have the n k ■ k — > singularity. We show in Appendix B.3.2 that the result takes 
the form 
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where /3y is the angle between the i-k and j-fc planes and the finite function F is given in 
Eq. (B.33) and is 0(l/t 2 ). 

5.2.3 Soft gluon inside measured jet k: Sfj Qas (T^) 

The contribution of a gluon emitted into jet 1 (the measured jet) from lines i and j is given 
by the integral 
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The singularity structure of this integral depends on whether or not k = i or j. Thus, we 
evaluate the case k = i or j and the case k ^ i,j separately below. 

Case 1: k = i or j We consider first <S™ cas (r*). Using the results of Sec. B.2 of the 
Appendix, we obtain the result in terms of tij, 
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Case 2: k ^ i,j The remaining contributions to the observed jet angularity are 5'™ cas 
for k i,j. Using the results from Sec. B.3.3 in the Appendix, this contribution is 
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(5.19) 



where G is 0(l/i 2 ) and is given in Eq. (B.36) and, again, Pij is the angle between the i-k 
and j-k planes. 
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Table 2: Summary of the divergent parts of the soft function at NLO. The first block contains the 
the observable-independent contributions: soft gluons emitted by jets i and j in any direction with 
energy E g < A in row 1; soft gluons entering jet k with E g > A (with k = i or j in the second row 
and k i,j in the third). In the last row of the first block, we summed over i and j and took the 
large- 1 limit to get the total 5 , " 1 ™ eas . Similarly, in the second block we give the contributions to 
the angularities (with k = i or j in the first row and k ^ i,j in the second) and summed over i 
and j and took the large-i limit to get S^ as in the third row. 

5.3 Total A^-Jet Soft Function in the large-i Limit 

In this section, we add together the necessary ingredients calculated above to obtain the 
iV-jet soft function from Eq. (5.8). The results for the divergent pieces are summarized in 
Table 2. In Sec. 6 we use Table 2 to show that the consistency of factorization is explicitly 
violated by terms of order 1/t 2 , and so in this section we give the full soft function (including 
the finite terms) to 0(l/t 2 ). 

Using color-conservation (^ i Tj = 0), we find that the observable-independent part, 
gunmeas ^ d en ; nec i [ n jt^ (5.9), is given for large t by 




(5.20) 
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We see that the finite parts of this contribution are sensitive to two scales, 2A and 2Atan R 



2 ■ 



For simplicity, in this paper, since we take tan(i?/2) ~ 0(1), we will choose only a single 
scale /ig to attempt to minimize logs in Eq. (5.20), where 

R 



Us = 2 A tan 



1/2 



2 ' 



(5.21) 



chosen as the geometric mean of the two. 

The remaining part of the soft function that is dependent on angularities as our choice 
of jet observable is the sum over 5*™® as (r*) (defined in Eq. (5.10)) for each jet i, where 
5™^ as (r*) is given by 
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(5.22) 



+ 0(l/t 2 ) . 



The finite part of this contribution is sensitive to the scale /i l s , where 



Ms 



OJiTX 



tan 1 ' - # 



(5.23) 



which, in principle, differs for each jet i and from the scale fig in the unmeasured part of 
the soft function Eq. (5.20). After discussing resummation of large logarithms through RG 
evolution, we will describe in Sec. 6.4 a framework to "refactorize" the soft function into 
pieces depending on multiple separated soft scales. This framework will enable us to resum 
logarithms of all of these potentially disparate scales. 



6. Resummation and Consistency Relations at NLL 

The factorized cross section Eq. (3.34) is written in terms of hard, jet, and soft functions 
evaluated at a factorization scale \i. Since the physical cross section is independent of /i, 
the anomalous dimensions of these functions are closely related. We derive and verify this 
relation in Sec. 6.3. In Sec. 6.1 and Sec. 6.2, we work out the form of the renormalization- 
group equations (RGEs) satisfied by the hard, jet, and soft functions, and obtain their 
solutions so that we can express each function at the scale [i in terms of its value at a scale 
//o where logarithms in it are minimized. In Sec. 6.4, we present a framework to refactorize 
the soft function and give the total resummed distribution in Sec. 6.5. 

6.1 General Form of Renormalization Group Equations and Solutions 

We will build solutions for the hard, jet, and soft functions from two forms of RGEs. The 
first form is for a function which does not depend on the observable r a and is multiplicatively 
renormalized, 



nbarc 



3.1) 
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and satisfies the RGE, 

M^(m)=7f(/x)F( M ), (6.2) 
where the anomalous dimension 7^ is found from the Z- factor by the relation 



and takes the general form, 



7f0) = - „ \ . li^-Z F (jj) , (6.3) 



V 2 

7f(m) = F F [a] In + j F [a] ■ (6.4) 



We call IV [a] the "cusp part" of the anomalous dimension and 7i?[a] the "non-cusp part". 
They have the perturbative expansions 



I >w = (£) I * + (£p* + -" («> 



2 



and 

'0!.,\ n . (a 



+ + (6 ' 6) 

The RGE Eq. (6.2) has the solution 

F(n) = U F (jji,lio)F( J i ), (6.7) 
where the evolution kernel LV is given by 

U F (»,, ) = e«^^y FM , (6.8) 

where -KV and ui F will be defined below in Eq. (6.15). 

The second form of RGE is for a function dependent on the jet shape r a and is renor- 
malized through a convolution, 

F barc (r a ) = J dr' a Z F (r a - <;/i)F(r» , (6.9) 

and satisfying the RGE 

/i~F(r a ; /x) = J dr' a lF (r a - r' a , fi)F(r' a ; /i) , (6.10) 
with an anomalous dimension calculated from 

/d 

dr' Z~ l (T a - T' a ; n) ^— Z F (T' a ; fj,) , (6.11) 

and taking the general form 

l F (T a ;n) = ~T F [a s ] (- \^] -In £s(T a )) + j F [a s ]S(r a ) . (6.12) 
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The solution of an RGE of the form Eq. (6.10) has the solution [82, 85, 86, 87, 52] 

F (i~a, M) = J dr> U F (r a - T' a ; fj,, no)F(T' a ; fJ-o) , (6.13) 
where the evolution kernel U F is given to all orders in a s by the expression 

UF{Ta > " ^ = Th^T w [ kF^J + ' (fU4) 

where j F is the Euler constant. 

In Eqs. (6.8) and (6.14), the exponents uj f and Kp are given in terms of the cusp and 
non-cusp parts of the anomalous dimensions by the expressions 
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(6.15a) 
(6.15b) 



In the case of Eq. (6.8) or if IV [a] happens to be zero, we define j F to be 1. To achieve NLL 
accuracy in the evolution kernels Up, we need the cusp part of the anomalous dimension 
to two loops and the non-cusp part to one loop, in which case the parameters uj f ,K f in 
Eq. (6.15) are given explicitly by 
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Here r 



Mm) 

Os(M0) 



and Po,/3i are the one-loop and two-loop coefficients of the beta function, 



da s 
djj, 
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(6.17) 



where (with Tr = 1/2) 



A-^-^ and ft^-HMi^C^. (6.18) 



The two- loop running coupling a s (fi) at any scale \x in terms of its value at another scale 
Q is given by 
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In Eq. (6.16), we have used the fact that, for the hard, jet, and soft functions for which 
we will solve, the cusp part of the anomalous dimension IV [a s ] is proportional to the cusp 



anomalous dimension r cusp [a s ], where 



r cusp [a s ] = (g) r° usp + (^) 2 r c 1 usp + • • • . (6.20) 

The ratio of the one-loop and two-loop coefficients of r cusp is [88] 

-^cusp ana ^ A are needed in the expressions of ujf and Kp for complete NLL resummation 
since we formally take a 2 . lnr a ~ 0(a s ). 

6.2 RG Evolution of Hard, Jet, and Soft Functions 
6.2.1 Hard Function 

The hard function is related to the matching coefficient Cat of the iV-jet operator in 
Eq. (3.33). If there are multiple operators with different color structures, Cm is a vec- 
tor of coefficients. The hard function is then a matrix H a f, = C$C a . The hard function is 
contracted in the cross section Eq. (3.34) with a matrix of soft functions. 

The anomalous dimensions of the matching coefficients C a have been calculated in the 
literature (for example, Table III of Ref. [89]). For an operator with N legs with color 
charges T 2 , the anomalous dimension is 

£ [T?r(«.) In £ + ^ K)| - \T(a s ) Ti-T; In ( "^ " *° + ) , 

i=l L 



(6.22) 

where Tj is a matrix of color charges in the space of operators, and 7« is given for quarks 
and gluons by 

3a s C F a s IIC A - 2N f 
lq = —ii , lg = R • 6-23 

The coefficient T(a s ) is the cusp anomalous dimension and is given to one-loop order by 
r(a s ) = a s /n. The anomalous dimension of the hard function itself is given by = 
Jq + 7Cjv) an d takes the form of Eq. (6.4), with cusp and non-cusp parts IV [a s ] and 
Jh [c*-s] given to one loop in Table 3 , with the result 

2 N 

lH (a s ) = -r(a s )T 2 ln^ - V 7i (a s ) - F(a s ) V T r T, In ^ , (6-24) 
where T 2 = Ei=i T ! is the 

sum of all the Casimirs, and the effective hard scale uh 
appearing as the scale oj in the logarithm in Eq. (6.4) is given by the color-weighted average 
of the jet energies, 
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Table 3: Anomalous dimensions for the jet and soft functions. We give the cusp and non-cusp 
parts of the anomalous dimensions, IV [a s ] and 7i?[a s ]. T is the cusp anomalous dimension itself, 
equal to a s /ir at one loop. 7, is given for quark and gluon jets in Eq. (6.23). The third column 
gives the value of jp appearing in Eq. (6.12) or Eq. (6.15). The last column gives the values of 
u> appearing in the logarithm ln/z 2 /u; 2 multiplying the cusp part of the anomalous dimension in 
Eqs. (6.4) and (6.12). The scale ujh is the color-weighted averages of all jet energies defined in 
Eq. (6.25). All rows except for the last are given in the large-i limit and in the last row we give 
the remaining terms that are present for arbitrary t. This last row explictly violates consistency at 
0(l/i 2 ). The first group of rows are needed for measured jets and the second group for unmeasured 
jets. In the large-t limit, for any number of measured and unmeasured jets, the consistency relation 
Eq. (6.33) is satisfied. 



6.2.2 Jet Functions 

There are two forms of jet functions, those for measured and those for unmeasured jets. 
Unmeasured jet functions J2' 9 (/u) satisfy multiplicative RGEs, with anomalous dimensions 
given by the infinite parts of Eqs. (4.17) and (4.26), 



7 j 1 =r(a s )T 2 ln 



cof tan 2 1= 



(6.26) 



which falls into the general form Eq. (6.4), with cusp and non-cusp parts of the anomalous 
dimension given in Table 3, and the scale to in Eq. (6.4) being Wjtan^. The part 7$ is 
given by Eq. (6.23). 

Measured jet functions satisfy RGEs of the form Eq. (6.10), with anomalous dimensions 
given by the infinite parts of Eqs. (4.11) and (4.25), 
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which takes the form Eq. (6.12) with cusp and non-cusp parts of the anomalous dimension 
split up as in Table 3, and the scale oj in Eq. (6.12) being Wj. 



6.2.3 Soft Function 

The total iV-jet soft function is given by Eq. (5.20) for unmeasured jets added to the 
sum over measured jets of Eq. (5.22). This soft function depends on the M jet shapes 
T a ) • • • j r Q Vf ' an d satisfies the RGE 

d f 
fi—Sin, t m ; fj)= / dr[ ■ ■ ■ dr' M 75 (n -t[,...,t m - r' M ; /x)5(r{, . . . , t' m ; y) . (6.28) 



d\x 

From the infinite parts of the soft function given in Table 2, we find the anomalous di- 
mension 7s (ti, . . . , tm', m) decomposes, as required by the consistency condition Eq. (6.33) 
given below, into a sum of terms, 

M 

7s(n, . . . , r M ; y) = 7 s nmcas (/^(n) • • • S( m ) + £ 7 s eas (^; m) II > ( 6 - 29 ) 

k=l j^k 

where the pieces 7^ nmeas (//) and 7™ eas (rfc; //) are given in terms of their cusp and non-cusp 
parts in Table 3, with the result 

7 — ( M ) = £ r(a s )T 2 In tan 2 - + T(a s ) £ T, • T, In ^ , (6.30) 

i v£j 

which takes the form of Eq. (6.4) with no cusp part, and 
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which takes the form of Eq. (6.12) with no non-cusp part, and the scale u in Eq. (6.12) 
being lo^j tan x ~ a 

The solution of the RGE Eq. (6.28) is given by 



. M 

S(n, . . . , r M ; 11) = / dr{ ■ ■ ■ dr' M S(r{, . . . , r' M ; /, )^ nmcas (/i, W>) II ^(^ " ^ M, A* 

fe=i 



(6.32) 



where [/^ nmeas is given by the form of Eq. (6.8) and Ug(r^) by the form of Eq. (6.14). 

The solution Eq. (6.32) is appropriate if all the scales appearing in the soft function 
are similar, and thus all large logarithms in the finite part can be minimized at a single 
scale /io- As we noted in Sec. 5.3, however, the potentially disparate scales WjT*, set by the 
jet shapes of the measured jets, and A, set by the cutoff on particles outside jets, appear 
together in the soft function, and logarithms of ratios of these scales may be large. In this 
case, the soft function should be "refactorized" into pieces depending only on one of these 
scales at a time. We describe a framework for doing so below in Sec. 6.4. 

But first, we verify the consistency of the anomalous dimensions for the hard, jet, and 
soft functions to the order we have calculated them. 
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6.3 Consistency Relation among Anomalous Dimensions 

We summarize the anomalous dimensions of the hard, jet, and soft functions in Table 3. We 
separate contributions to the jet and soft anomalous dimensions that arise from measured 
jets, from unmeasured jets, and those that are universally present. In all rows except the 
last row, we take the large-i limit and give the additional terms that arise (from the soft 
function) for arbitrary t. 

Consistency of the effective theory requires that the anomalous dimensions satisfy 

o = (7*00 + 7 r mcas (^) + E 7j,(m))^) + E (j^i; m)) + 7r as (r: ; . 

i^meas «6mcas 

(6.33) 

From the results tabulated in Table 3, up to corrections of 0(l/t 2 ), we see that, remark- 
ably, this relation is indeed satisfied! This is highly nontrivial, as jet and soft anomalous 
dimensions depend on the jet radius R and the jet shape r a , while the hard function does 
not; in addition, the hard and soft functions know the directions rii of all iV jets, while the 
jet functions do not. These dependencies cancel precisely between the ^-dependent pieces 
of unmeasured jet contributions to the jet and soft functions, between r a -dependent pieces 
of the measured jet contributions, and between m ■ n^-dependent pieces of the hard and 
soft functions. The sum of all jet and soft anomalous dimensions then precisely matches 
the hard anomalous dimensions, satisfying Eq. (6.33). 

Making the satisfaction of consistency even more nontrivial, individual contributions 
to the infinite part of the soft function, and therefore its anomalous dimension, given by 
Table 2 depend on the energy cut parameter A as well. However, these terms cancel in the 
sum over the contributions Sfj and S\j in the first two rows of Table 2. The double poles 
of the form Mn A arise from regions of phase space where a soft gluon can become both 
collinear to a jet direction (giving a 1/e) and soft (giving a In A). These regions exist in 
the integral over all directions giving S 1 ^ but are subtracted back out in the contributions 
In the surviving "Swiss cheese" region the regions giving these double poles are cut 

out. 

The 0(l/t 2 ) terms that violate consistency arise whenever there are unmeasured jets. 
While this limit is not required for the contribution of measured jets to the anomalous 
dimension to satisfy the consistency condition Eq. (6.33), the finite parts of measured 
jet contributions to the soft function contain large logarithms of ui/A that can not be 
minimized with a scale choice but are suppressed by 0(l/t 2 ) (cf. Eq. (B.37) of Appendix 
B). This is the manifestation of the fact that jets need to be well-separated, as explained 
in Sec. 3. For the remainder of the paper, we work strictly in the large-i limit. 

It may seem mysterious that the calculations of the hard, jet, and soft functions them- 
selves and requiring their consistency lead to the condition of a large separation parameter 
t. Although we already specified qualitatively in the proof of factorization the requirement 
of well-separated jets, it may not be immediately apparent where it is implemented in the 
actual calculations. It enters in the definition of the collinear jet functions. In the large-f 
limit, the N jets are infinitely separated from one another according to the measure given 
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by Eq. (1.4). And indeed, when iV-jet operators are constructed in SCET, each collinear 
jet field contains a Wilson line W n , defined below in Eq. (3.9), of collinear gluons in the 
direction n that were emitted from the back-to-back direction n, for which the separation 
measure t — > oo. (This is similar to QCD factorization proofs of hard scattering cross 
sections, e.g. in [17], in which this direction n is chosen to be along some arbitrary path 
£ that is separated by an 0(1) amount from the jet direction n.) Furthermore, the rij- 
collinear jet function knows only its own color representation, and not those of the other 
jets. Meanwhile, the hard and soft functions we calculate "know" about all N jets and 
their precise directions and color representations. Therefore it is no surprise that, when 
we actually calculate the anomalous dimensions of these functions, we find that they are 
consistent with one another only in the limit that the separation parameter t — > oo. 

6.4 Refactorization of the Soft Function 

Our results for the soft function in Sec. 5.3 make clear that in general there are multiple 
scales that appear in the soft function: the fig, . . . ,/iff associated with the M measured 
jets and the scale jKo associated with the out-of-jet cutoff A (see Eq. (5.21)). When these 
scales are all comparable, we can RG evolve the soft function from the single scale fis- 
However, when any of them differ considerably from the others, we need to "refactorize" 
the soft function into multiple contributions, each of which is sensitive to a single energy 
scale. For illustration, take the scales fig to be such that fig ^> fig 3> • • • 3> /iff as in 
Fig. 7. We also take fJ'g 1 <C fig <C fig for our discussion, but the result is independent of 
whether fig lies in the range fi s <fig< fig 1 or not. 

We can express the soft function appearing in Eq. (3.34) as 



S{r l a ,rl r a M ; fi) = (O|0<jG(A - A) JJ 5(r* - f<)O s |0> , (6.34) 



Ms~ 

m|- 



Ms- 
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where the operator r* returns the contribution to r a of final-state soft 
particles entering jet i, and A returns the energy of final-state soft particles 
outside of all TV" jets. We have kept the dependence on the scales fi l s and 
on A implicit on the left-hand side. Figure 7: Soft 

There are N Wilson lines appearing in the operator Os, scales. 

O s = Y 1 ...Y M Y M+l ...Y N , (6.35) 

corresponding to the M measured jets and N — M unmeasured jets. The scales associated 
with soft gluons entering the M measured jets whose shapes are measured to be n, . . . , tm 
are fig, . . . , fig 1 , given by Eq. (5.23). The scale associated with soft gluons outside of 
measured jets is fig given by Eq. (5.21). We have illustrated the ladder of these scales 
in Fig. 7. Each of these soft scales can be associated with different soft fields A\ whose 
momenta scale as Xf uji where Aj is associated with the typical transverse momentum AjWj 
of the collinear modes for the ith jet. For measured jets, Aj is determined by r*, while for 
unmeasured jets Aj ~ tan(i?/2). For soft gluons outside jets, however, the soft momentum 
is set by the cutoff scale A, which is why fig appears in the ladder of Fig. 7. 
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At a scale \i larger than all /j, s and fig, the soft function is calculated as we presented 
in Sec. 5. In particular, we take the quantities r* and A to be nonzero and finite. At a scale 
\i below /ig, we integrate out soft gluons of virtuality fi s and match onto a theory with soft 
gluons of virtuality fj, s . The scale n s associated with is taken to infinity, and phase space 
integrals for soft gluons entering the measured jet 1 become zero (see, e.g., Eq. (B.17)). 
Therefore, the matching coefficient from the theory above /j, s to the theory below is just 
the measured jet 1 contribution S' meas (r^) to the soft function given by Eq. (5.22). The 
same occurs when matching from the theory above each scale /j, s for soft gluons entering 
measured jet i to the scale below /j, s , giving a matching coefficient S' meas (r^). 

When going through the scale fj,g, in the theory above this scale, the calculation of 
unmeasured contributions to the soft function gives the result Eq. (5.20), by treating A as a 
nonzero, finite cutoff. In the theory below fig, we take A to infinity, making all phase space 
integrals originally cutoff by A to be scaleless and thus zero. So the matching coefficient 
between the theory above and below fig is just < 5 mnmeas . 

After performing the above matchings all the way down to the lowest soft scale in 
Fig. 7, we find that the original soft function S(r^, . . . , t^ 1 ; /x) can be expressed to all 
orders as 

M 

S(r a \ . . • , r a M ; //) = S" (/x) ]J S mcas (ri; /x) <0| O^Os |0> , (6.36) 

i=l 

where to next-to-leading order S me£LS and g unmeas are given by 

S unmcas (fi) = 1 + Sfl) mcas {fi) (6.37) 
5 meas (r^ /x) = 5(r*) + S^"(r*; /i) , (6.38) 

where S^ cas is given by Eq. (5.20) and 5^ as is given by Eq. (5.22). Note that no 
operators restricting the jet shape or the phase space appear in the final matrix element of 
soft fields living at the lowest scale on the ladder in Fig. 7. If all the scales on the ladder 
are at a perturbative scale, we can now just use (O s Os) = 1 to eliminate the final matrix 
element. If any scale is nonperturbative, we should have stopped the matching procedure 
before that scale, and defined the surviving soft matrix element still containing additional 
delta function operators as a nonperturbative shape function. 

Since the factors S unmeas (n) and 5' meas (r*, fi) are now matching coefficients between 
two theories above and below the respective scales fig and /j, s , we can run each of the 
individual factors separately from their natural scale, instead of from a single soft scale 
as in Eq. (6.32). The result for the RG-evolved soft function is then Eq. (6.36) where each 
factor at NLO is given by the solution of its RGE, 

^unmeas^ = f/unmcas^ ^A^unmcas^A) 39^ 

5 meas (r:, M ) = J dT'U s K - r'; /x, // 5 )S mcas (r', /£) . (6.39b) 

These solutions allow us now to resum logarithms of all of the scales appearing in the 
ladder in Fig. 7 when these scales are widely disparate. However, the result we obtained 
in Eq. (6.28) when we took all scales to be of the same order and had a single soft scale 
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has the form Eq. (6.39) at NLL accuracy. We will use equation Eq. (6.39) in all cases to 
interpolate between these two extremes. 

6.5 Total Resummed Distribution 

Collecting together the above results for the running of hard, jet, and soft functions in 
the factorized cross section Eq. (3.34), we obtain the RG-improved iV-jet cross section 
differential in M jet shapes, 
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where ojh is defined by Eq. (6.25), the evolution parameters wf(w, /j-f) an d Kf{^-,^f) are 
defined in Eq. (6.15), and we have defined the collective parameters, 
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Using results from Appendix C, we obtain the functions /} s generated by the finite pieces 
of the measured jet and soft functions, 
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where Cj = 3/2 for quark jets and /?o/(2Ca) for gluon jets. r™ ax is the upper limit on r* 
found in the finite part of the naive jet function, given in Appendix A. H(—l — Qi) is the 
harmonic number function, with fij given by Eq. (6.41b). tp^ 1 ' is the first derivative of the 
digamma function, ^ l \z) = (d/dz)[T'(z)/T(z)}. The terms d, ls are additional contribu- 
tions from the finite parts of jet and soft functions that do not contain any logarithms, 
where dg = — 7r 2 /24, and dj is given in Eq. (C.6) in the Appendix, djg are not needed 
at NLL accuracy. Similarly, the terms containing large logarithms in the unmeasured jet 
functions and unmeasured contribution to the soft function are 
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where dj is the part of the unmeasured jet function containing no large logarithms (given 
in Eqs. (A. 19) and (A. 30) in the Appendix). 

The finite parts of the measured and unmeasured jet and soft functions given in 
Eqs. (6.42) and (6.43) show that to minimize large logarithms in the 0(a s ) finite parts in 
the resummed distribution Eq. (6.40), we should choose initial scales for the running to be 
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These choices eliminate all large logarithms in the 0(a s ) hard, jet, and soft functions. 
They still leave logs of tan ^ and n, • rij in the unmeasured part of the soft function, and 
logs of tan 7^ in the measured jet function, but we already take R numerically of 0(1) 14 
to minimize power corrections from our implementation of the jet algorithm as discussed 



in Sec. 3.4, and rij 



1 since the jet separation parameter tij is large compared to 1. 



All logs of R, A, and r* are eliminated in the unmeasured jet function and measured part 
of the soft function. 



14 We still consider tan(7?/2) to be of order Afe in the collinear sectors describing unmeasured jets, as 
implied by Eq. (6.44). This means Afe is effectively much larger than the parameter Ai in a measured 
jet sector. In fact, note that Eq. (6.44) tells us that tan ^ must be parametrically larger than (r„) 2 ~ a ; 
otherwise, the jet scale falls below the soft scale in the measured jet sectors, invalidating the use of SCET 
and, thus, the validity of the factorization theorem. 
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7. Plots of Distributions and Comparisons to Monte Carlo 



Having resummed the jet shape distributions in r a to NLL accuracy, in this section we 
plot the distributions for various values of a and R, compare to Monte Carlo simulated 
events, and perform scale variation on the resummed distribution. We use the process 
e + e~ — > 3 jets to study our predictions of jet shapes, where the jets arise from partons 
in the "Mercedes-Benz" configuration, with each parton having equal energy. In these 
configurations, the partons lie in a plane and are equally separated with a pairwise angle 
of 27r/3. This allows us to study event shape distributions of well-separated jets where t is 
reasonably large. We choose three values of R to study, R = 1.0, 0.7, and 0.4. With these 
values of R, the 1/t 2 suppression factor for corrections to the large-i limit are 0.10, 0.044, 
and 0.014 respectively. We will measure the angularity of only one of the three jets; the 
other two jets will be unmeasured. 

In general, the T, • Tj color correlations in the soft and hard functions lead to operator 
mixing in color space under RG evolution. This implies that the RG kernels Us and Uh 
are matrices in color space and must be studied on a process- by-process basis (see, e.g., 
[89, 90, 91, 92, 93, 94]). For the case of N = 2, 3 jets there is only one color-singlet operator 
and hence no mixing. This can be seen, for example, by noting that all color correlations 
reduce to the Casimir invariants (Cf and Ca) in this case (cf. Appendix D). We have 
restricted the example process we use in this work to N = 3 jets, avoiding the additional 
complication of color-correlations that comes with a larger number of jets. 

The NLL resummed distribution for one quark or gluon jet shape (jet 1) in a three-jet 
final state, written as the derivative of the radiator Eq. (1.5), is 
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where the various evolution parameters iOj s , fC are all defined in Eqs. (6.15) and (6.41), 
and fj t s are given by fj t s in Eq. (6.42) with the dj t s terms set to zero (accurate to NLL). 
The best scale choices Eq. (6.44) for this case are 
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In Eq. (7.1) we have used tree-level initial conditions for the hard, jet, and soft functions 
at these scales. Eq. (7.1) evolves these functions to the arbitrary scale // at NLL accuracy. 
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Figure 8: Quark and gluon jet shapes for several values of a and R. The NLL resummed distribu- 
tion in Eq. (7.1) is plotted for a = — |, — |, 0, |, ^ for quark and gluon jets with R = 1, 0.7, 0.4. The 
plots are for jets produced in e + e~ annihilation at center-of-mass energy Q = 500 GeV, with three 
jets produced in a Mercedes-Benz configuration with equal energies Ej — 150 GeV, and minimum 
threshold A = 15 GeV to produce a jet. 



With these choices, we plot Eq. (7.1) in Fig. 8 for several values of a and R for a quark 
or gluon jet shape in a three-jet final state in e + e~ annihilation at center-of-mass energy 
Q = 500 GeV. 15 The jets are chosen to be in a Mercedes-Benz configuration, where all 
jets have equal energies of 150 GeV. We choose the jet energy cutoff A to be 15 GeV. We 
choose the factorization scale to be [i = fin- 

15 The distributions plotted with the /j,s terms included in Eq. (7.1) exhibit a small negative dip near r a = 
(not shown) that can be cured by convolving with a nonperturbative shape function with a renormalon- 
free gap parameter [38, 54]. This is beyond the scope of the present work, so we only plot the perturbative 
distributions where they are positive. 
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Figure 9: Quark vs. gluon jet shapes with comparison to Monte Carlo. Solid, straight curves 
represent the resummed jet shape distribution in Eq. (7.1), and jagged curves are histograms from 
the Monte Carlo, normalized as described in the text. The solid histogram has no hadronization, 
while the dashed histogram includes the effects of hadronization. The distributions are plotted 
for a = — ^,0, \ with quark (blue) and gluon (red) jets compared on the same plot, for jets of 
size R — 1.0,0.7,0.4. Gluon jet shape distributions peak at larger values of r a than quark jets, 
indicative of their broader shape. The plots are for jets in e + e _ annihilation at center-of-mass energy 
Q = 500 GeV, with three jets produced with equal energies Ej = 150 GeV, angular separation 
i\) ~ 2tt/3 between all pairs of jets, and minimum threshold A = 15 GeV to produce a jet. 



We compare the results of a jet algorithm implemented on Monte Carlo simulated 
events with our NLL resummed predictions for a variety of a and R values in Fig. 9. Because 
the resummed NLL distribution we choose to study applies to an exclusive process, three- 
jet events in the Mercedes-Benz configuration, we implement cuts on the simulated events 
to obtain a sample that matches onto this configuration. We use MadGraph/MadEvent 
v. 4. 4. 21 [95] to generate parton-level e + e~ — > qqg events at a center-of-mass energy Q = 
500 GeV, with cuts imposed to obtain partons in the Mercedes-Benz configuration. We 
shower and hadronize the events with Pythia v. 6. 414 [96] using p^-ordered parton showers. 
The process of hadronization will induce a shift in the angularity distribution, which we do 
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not model in our resummed distribution. Therefore, we produce two samples: one sample 
with only QCD final-state showering, no initial-state radiation, and no hadronization, and 
another sample with complete showering, initial-state radiation, and hadronization. The 
anti-kx jet algorithm is run on the final state particles from Pythia, and we use Fast Jet [97] 
to implement the jet algorithm. The anti-kT algorithm is particularly well suited for this 
comparison, as very few particles at an angle 8 > R to the jet axis are included in the jet. 
With anti-kx, the phase space cut on an individual particle matches well with the phase 
space cuts in the next-to-leading order calculation. 

To select a sample of events to compare to our resummed distributions, we make cuts on 
the final state jets, requiring each of the three hard, well-separated partons from MadGraph 
to be associated with a jet. This involves a cut on the jet direction and momentum: 

Pparton • Pjct > Q Q ^ | |p P arton| ~ |Pjet| | < q j Ej (7.3) 
|Pparton| |Pjct| | Pparton | 

We analyze events passing these cuts, and tag each associated jet as coming from a quark 
or a gluon based on which parton it matches onto. The angularity value for each jet is 
computed from the constituent particles of the jet, using the matching parton direction as 
the jet axis. The jet direction only differs from the parton direction by a power correction 
(see Sec. 3.2). In Fig. 9, we isolate some of the quark and gluon jet shapes in Fig. 8 and 
compare to Monte Carlo events. 

The relative normalization between the distribution of Monte Carlo events and the 
NLL resummed angularity distribution requires some explanation. Both our calculation 
and the Monte Carlo simulation are most accurate in the region that includes the peak 
of the distribution and the larger-r side of the peak, but both are inaccurate as r — > 
and in the tail region. Therefore, each will differ in the relative normalization between 
the peak region and the tail region. An accurate prediction of the tail region requires 
matching onto a calculation at fixed-order in a s in full QCD as in [43, 53, 54]. In Fig. 9, we 
choose to normalize the area of the Monte Carlo distribution to the total area of the NLL 
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Figure 10: Location of peak of jet shape distribution as a function of a for quark and gluon jets. 
We plot the value of r a at the peak of the jet shape distribution for a between -1.0 and 0.8 for 
quark vs. gluon jets, using R = l, 0.7, 0.4. 
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Figure 11: Scale variation of quark and gluon jet shapes. For a = and R = 0.7, we display 
the variation of the NLL resummed jet shape distributions with the hard scale fin, the jet cutoff 
scale //g, the unmeasured jet scales fij , the measured jet scale fjtj(r a ), and the measured soft scale 
Hs(T~a)- In each case we vary the scale between 1/2 and 2 times the natural choices in Eq. (6.44), 
except for the measured soft scale, which we varied between 1 and 2 times the choice in Eq. (6.44). 
We keep the factorization scale fixed at the default hard scale given by Eq. (7.2), \i = lu.- l . 
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resummed theory distribution. We find the area under the theory curves for quark and 
gluon jets to be approximately 0.3 for R = 0.4, 0.5 for R = 0.7, and 0.7 for R = 1. A more 
accurate prediction of the normalizations may require summing remaining unsummed logs 
of the phase space cuts in the theory and Monte Carlo predictions. These plots should be 
interpreted as comparisons of the predictions of the shapes in r a and these shapes' scaling 
as we vary a and R, rather than the overall normalization. 

The shapes of the theory and Monte Carlo distributions are largely similar, though they 
display noticeable differences at the leftmost endpoint near r a = and in the "sharpness" of 
the peak. These may be due to the different ways the two approaches deal with the growth 
of the strong coupling for small r a , the different orders of log resummation (LL vs. NLL) 
and the need to match the tails onto fixed-order QCD predictions. Since neither the Monte 
Carlo nor theory partonic predictions without inclusion of hadronization effects is yet a 
prediction of a physically observable quantity, we use this comparison as an intermediate 
diagnostic rather than a conclusive test of either method. Nevertheless, comparing the way 
the shapes of the distributions and locations of the peaks vary over the range of values of 
a and R that we sample, the behavior agrees very well between the theory distributions 
and the Monte Carlo distributions without hadronization for both quark and gluon jets. 

In Fig. 10 we plot the location of the peak of the jet shape distributions as a function 
of a for three values of R, displaying the different variation of the peak of quark and 
gluon jet shape distributions. The peak value increases with increasing R and a, as wide 
angle radiation is included (increasing R) and less suppressed (increasing a). Although 
the difference in the peak value between the quark and gluon jet angularity distributions is 
large, the width of each distribution creates substantial overlap in angularity values between 
quark and gluon jets. Distinguishing between quark and gluon jets using jet angularities is 
a complex task which we will explore in future work; for now, we note only that the NLL 
resummed distributions indicate that discrimination between quark and gluon jets using 
jet angularities is possible. 

As a rough estimate of the theoretical uncertainty in our NLL resummed predictions, 
we show in Fig. 11 the change in the a = quark and gluon r a jet shape distributions for 
R = 0.7 when the various scales that appear in the resummed cross section Eq. (7.1) are 
varied. These are the initial scales at which the hard, jet, and soft functions are evaluated 
to minimize logarithms in the NLO fixed-order part, from which the evolution kernels run 
them to the common factorization scale [i. In the top row of Fig. 11, we vary \x between 
u>#/2 and 2u>h- The tiny variation is a sign of the consistency condition satisfied by the 
anomalous dimensions in Eq. (6.33). In the next four rows, we vary the hard scale /U#, 
the soft jet energy cutoff scale fig, the unmeasured jet scales /ij 3 , and the measured jet 
scale /ij(r^) between half and twice the natural values given in Eq. (7.2). In the last row, 
we vary the measured soft scale /^(t*) between one and two times the value in Eq. (7.2). 
This is because too low a value of /^(t*) as r a — > brings it into the nonperturbative 
region where a s (fig) blows up, so that the perturbative estimate of uncertainty is not so 
meaningful. We note that, while the uncertainty in the vertical scale of the distributions 
is considerable in some cases, the location of the peak is much more stable. 

Finally, in Fig. 12 we give a sense for how robust our theoretical predictions are for other 
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Figure 12: Jet shapes for other kinematic configurations. We compare our theoretical predictions 
to Monte Carlo simulations for the shape Tg (a = 0) for a quark or gluon jet found in a three-jet 
configuration where the two jets with narrowest separation angle V'ncar have equal energy. We 
consider the two cases linear = 7r/2 and tt/3. In the first row, we plot shapes of one of the jets in 
the "near" pair. The blue solid curve is the shape of quark jet found near a gluon jet, the green 
dotted curve is a quark found near an antiquark, and the red solid curve is a gluon found near a 
quark. In the second row, we compare shapes of a quark or gluon jet found far from the near pair. 



kinematic configurations. We consider e + e~ — > qqg events where the angle V'near between 
two partons is either it/ 2 or tt/3, and these partons have equal energy. We find jets using 
the anti-kT algorithm with R = 0.4, and plot jet shapes for a = 0. The selection cuts to 
choose events from the Monte Carlo are the same as the Mercedes-Benz configuration. In 
these events there are five distinct characterizations for a single parton. If the event has 
the quark (or antiquark) as the "far" (most well separated) parton, then each parton in 
the event is distinct: there is the far quark, the near quark, and the near gluon. If the 
event has the gluon as the far parton, then there are only two distinct types of partons: 
the far gluon and the near quark (antiquark). In Fig. 12, we plot all these configurations 
for both Vncar = tt/2 and V'near = tt/3. The agreement between the theory predictions 
and the Monte Carlo are as good as in the Mercedes-Benz case, a good indication that 
our calculation applies to a broad range of kinematic configurations of multijet events. 
Additionally, we observe features consistent with our intuition about the relative differences 
between the jet shape distributions between different jets in these configurations. As one 
would expect, the distribution of near jet shapes is weighted more heavily towards larger 
r a than the far jet shapes, due to the enhanced soft radiation in the near jet system. When 
the near quark is near a gluon, the distribution is weighted more heavily towards larger 
r a than when the near quark is near an antiquark, due to the enhanced radiation coming 
from a gluon rather than a quark. These distributions serve as further evidence that jet 
shapes may be an effective discriminant between quark and gluon jets. 
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8. Conclusions 



In this work, we have factorized an iV-jet exclusive cross section differential in M < N 
jet observables and resummed global logarithms of the jet observable r a to NLL accuracy, 
leaving summation of non-global logarithms to future work. We demonstrated that the 
anomalous dimensions of the hard, jet, and soft functions in the factorization theorem 
satisfy the nontrivial consistency condition Eq. (6.33) to 0(a s ), for any number of quark 
and gluon jets, any number of jets whose shapes are measured, and any size R of the jets, 
as long as the jets are well-separated, meaning t S> 1. This condition ensures the validity 
of an effective theory with N collinear directions that are assumed to be distinct. We 
identified and estimated important power corrections to the factorized form of the cross 
section. We also illustrated that zero-bin subtractions give nonzero contributions to the 
anomalous dimensions crucial for consistency. 

Armed with consistent factorization and the fixed-order jet and soft functions, we 
resummed large logarithms in the jet shape distribution by running each individual function 
from the scale where logs in it are minimized to the common factorization scale fi. We 
thereby resummed, to NLL accuracy, global logs of the jet shape r a and logs of the scale 
A/Ej of soft radiation outside of jets, but leaving some non-global logs and logs of the 
angular cut R (but we took R to be numerically of order 1) . This is the first such calculation 
of a resummed jet shape distribution in an exclusive multijet cross section. 

We constructed a framework to deal with all the scales that appear in the multijet 
soft function which depends on the values r* of all M jet shapes and the phase space cuts 
A, R. By refactorizing the full soft function into individual pieces depending on one of 
these scales at a time, we were able to sum logs of ratios of these scales. 

We demonstrated the accuracy of our results by comparing our resummed prediction 
for quark and gluon jet shapes in e + e~ — > 3 jets to the output of Monte Carlo event 
generators, MadGraph/MadEvent and Pythia. We compared our predictions with the 
Monte Carlo output without hadronization. The changes in shape and location of the peak 
value as functions of a and R match quite well between the theory and Monte Carlo. 

Our results provide a basis for future studies of other jet observables at both e + e~ and 
hadron colliders, requiring recalculation of those parts of our jet and soft functions that 
depend on the choice of observable. Studying jets at hadron colliders requires constructing 
observables appropriate for that environment and the switching of two of our outgoing jets 
to incoming beams, which can be described by beam functions in SCET [62]. 

Precision calculations of jet shapes will allow improved discrimination of jets of different 
origins. We are applying the results of our predictions of light quark and gluon jet shapes 
to distinguish quark and gluon jets with greater efficiency than achieved before. Extensions 
to the shapes of heavy jets and calculations of other types of jet shapes such as the ^(r/R) 
shape introduced in [14, 15, 16] can also be performed. 

Note added in final preparation: As this paper was being completed, Ref. [98] appeared 
reporting the calculation of a quark jet function for a jet defined with a Sterman- Weinberg 
algorithm and whose invariant mass s is measured. This jet function is related to our 
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measured jet function J2(r a ) for a cone jet at a = given in Eq. (4.11), since s = u 2 tq. 
We have checked that the corresponding results between the two papers agree. 
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A. Jet Function Calculations 

A.l Finite Pieces of the Quark Jet Function 

Measured Quark Jet Function The finite pieces the jet functions, which depend on 
the jet algorithm, share common features. For cone-type algorithms, the finite piece of the 

J q 



naive part of the quark jet function, Jli„(T a ), is given by 
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where in this Appendix, plus distributions are defined by [62] 



[e(x)g{x)]+ = lim— [Q(x - e)G(x)], with G(x) = / dx'g(x'), (A.2) 



defined so as to satisfy the boundary condition J dx[Q(x)g(x)]+ = 0. The quantity Z< 
depends implicitly on r a and R and is given by 
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Jx conc x x conc z 

The parameter x con e = ^cone(^a) is the lower limit on the x = q~ /oj scaled gluon momentum 
integral imposed by the cone restriction. It is given by the solution of the equation 
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where / CO ne(^) is defined as 

/conc(x) = X 2 ~ a [x- 1+a + (1 - X y 1+a ] (A.5) 

in the range < x < 1/2, which is plotted in Fig. 13A. The limit r™ ax is given by the 

maximum value over x of Eq. (A.5). Similarly, for kT-type algorithms, Jj? (r a ) is given by 
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Figure 13: Regions of integration for the (A) cone and kr-type algorithms for (B) a > — 1 and 
(C) a < —1. The allowed region of x is when the (blue) functions / CO ne, k T ( x ) ue above the (red) 
lines of constant T a j tan( 2 ~ a ) R/2. When a < — 1 for the kx algorithm, there are two regions of 
integration when r a > 2 a ~ 2 tan( 2 ~ a ) R/2. 
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is satisfied. We plot this region in Fig. 13B and C for the cases a > — 1 and a < — 1, 
repsectively. The boundaries of this region are the points X\$, illustrated in the figure, and 
are given by the equation 

A T (xi, 2 ) = tan £^, (A.9) 

where we take X2 > x\ if X2 exists. The upper limit r™ ax is given by the maximum value 
over x of the right-hand side of Eq. (A. 8). In general, the constraint Eq. (A. 8) is symmetric 
about x = |, and so the region 1Z is symmetric about the same point. In general, if a > — 1 
or T a < 2 a ~ 2 tan( 2-a ) then 1Z is a single range in x. Otherwise, 1Z is two disjoint ranges 
in x. Since r a > 2 a_2 tan( 2_a ) ^ can ° m y occur for a < — 1, we can write 2^ as 
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Note that Z?onc and Zj? involve the same integrand, but for each algorithm the integral 
is over different ranges. In addition, both x conc and x\ approach the same limiting value 
for small r a , 

(All) 



X > 



tan( 2 " a ) f ' 



Thus, we can extract the small r a behavior of both distributions by writing 



r a \ x 



In 



1 — x 



tan( 2 " a ) f x 



In 



tan( 2 " a ) f 



(A.12) 
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where x = x cone or x\ for the cone and kx algorithms, respectively. Defining 

1 — x 

r q (x) = 3x + 2 In , 



(A.13) 



using Eq. (A. 12), and including the zero-bin subtraction in Eq. (4.10), we find that the 
finite distributions of the full measured quark jet functions are 
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and 
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For a = 0, these expressions for the jet functions can be simplified further to give 
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for the cone jet function, and 
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for the kx jet function. In Eq. (A. 15b), x\ is given by its value for a = 0, 
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In Eq. (A. 15), we have divided the cone and kp jet functions into the contribution </^ c j(to) 
to the inclusive jet function [75, 76], given by 



2 oj 2 uj 2 2 2 



e(r ) 



TO 



3 , a 2 
- + 2 In -^r 
2 w 2 - 



(A.17) 

and algorithm-dependent parts. The algorithm-dependent part of the a = cone jet 
function Eq. (A. 15a) agrees with [98]. Note that if one takes R to be parametrically larger 
than To (cf. Sec. 3.4 and Eq. (6.44)), the algorithm-dependent parts of Eq. (A. 15) are 
power suppressed, and the cone and kp jet functions reduce to the inclusive jet function. 

Unmeasured Quark Jet Function The finite pieces for the unmeasured quark jet 
function are 
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where the constant term cij alg is given by 
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for the cone and kp algorithms, respectively. 

A. 2 Finite Pieces of the Gluon Jet Function 

Measured Gluon Jet Function The finite distributions of the naive gluon jet function 
are given by 



Jcone(Ta) = S(T a ) 
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and 



JL(r a ) = S(r a 
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where the integrals If, are given by 
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(A.21) 
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with the cone and kp regions of integration the same as for the quark jet functions. The 
value r™ ax is the same as in the measured quark jet function, for the respective jet algo- 
rithm. 
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Going through similar steps as for the quark jet function, defining 



r g {x)=2C A ln 
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(A.23) 



and using Eq. (A. 12) to make all logarithmic dependence on r a explicit, we find for the 
cone and kp-type jet function finite distributions 



J conc( r a) = S(r a ) 
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where x cone and 3^2 are given in Eqs. (A. 3) and (A. 8). 

For a = 0, the simplified result for the gluon cone jet function is 



^coneOo) = ^fncl( T o) + 



G(r )9 (tan 



2 i? 



ro 



r + tan 2 f 



r + tan 2 f 



+ 



9 (r - tan 2 f ) 



T() 



(A.25) 



where 



f{x) = C A i-x A - x + 4 - T fl JV) -x^ - 2x + 2 



(A.26) 
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For a = 0, the gluon k/p jet function is given by, 

G(T )eQt a n 2 f -r ) 
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where cci is given by Eq. (A. 16), and r g {x\) is given by Eq. (A. 23). In Eqs. (A. 25) and 
(A.27), the contribution Jf nc ]{jo) to the inclusive gluon jet function [34, 35, 78, 79] is 
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As for the quark jet functions Eq. (A. 15), the gluon jet functions split up into the inclusive 
jet function and algorithm-dependent pieces that are power suppressed for tq <C R. 

Unmeasured Gluon Jet Function For the unmeasured gluon jet functions, the finite 
pieces are given by 
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where the constant part d 9 j alg for the cone and kT algorithms is given by, respectively 
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B. Soft function calculations 



B.l Sjf 1 



To evaluate the expression Eq. (5.11), we first define 

\ - e r(l - e) 2vr V 4A 2 ) ±l l]± ^ ™ j) ' 



(B.l) 



We need X incl to O(e). Working in a coordinate system with Hi aligned along the z-axis 
and Hj in the xz- plane and defining n = 1 — 7ij • rij = n z - , we have 
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where z = ^~7f~^p^ " • The integration over u = cos 6 has singularities at the points u = 1 
and u = n which correspond to z = 1 and 2 = 0, respectively. To isolate these singularities, 
we split the integration over u into the ranges (—1,5) and (5, 1) where n < S < 1, 



-*i( n . . % .) = Zf cl (ni • nj ) + X^\m ■ rij) , (B.3) 



where 
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Over the range of integration of u in Z{ ncl , z G [0, 1) for 5 < 1. For X™ 01 , 2 £ (0, 1]. 

Furthermore, the singularity at u = n in lj ncl is made more explicit through the use 
of the identity 



2 F 1 [ 1,1-1- e;z) =f a {z) + h{z) 
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/ a (z) gives an 0(1/ e) contribution and we proceed by using the following trick that we 
exploit multiple times throughout the Appendix. 

To integrate a product of functions f(x,e)g(x,e) where / diverges at the point xq as 
(x — xq)~ 1+ci ( € \ we write the integation as 

dx f(x, e)g(x, e) = j dx f(x, e)g(x , e) + J dx f(x, e) (g(x, e) - g(x , e)) . (B.6) 

The first integral has relatively simple x dependence since g(xo,e) does not depend on x. 
The term in parenthesis in the second integral vanishes as x — xq for regular functions g 
and so the entire integrand can be expanded in e. 

We can now evaluate f a (z) by adding and subtracting the non-singular part of the 
integrand (which is the hypergeometric function) evaluated at u = n as in Eq. (B.6), 
whereas fb(z) is 0(e) and so we can simply expand about e = 0. Adding these contributions, 
we find that 
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For 2, ncl , the part of the integrand that is not singular at u = 1 is everything that 
multiplies (1 — and so we add and subtract this part as in Eq. (B.6). This gives 
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The integrals in Eqs. (B.7) and (B.8) give rise to many terms. However, we find that, 
after some lengthy algebra, the dependence on <5 cancels in the sum as it must and that 
the result can be simplified to 
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B.2 and S™^) 
B.2.1 Common Integrals 

In evaluating the soft contributions S\a and S^ eas (T^), we find an integral of the following 
form: 
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To evaluate this integral, we add and subtract the part of the integrand that is not singular 
at u = 0, namely f(u; f3,t), as in Eq. (B.6). This allows us to write 
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and that the expansion of the hypergeometric about e = for t > u is 
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Evaluating the integals, we obtain 
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B.2.2 S^(ri) 

To evaluate Eq. (5.17) for the case that k = i, we use light cone coordinates in the frame 
of jet i, k + = rii ■ k and k~ = hi ■ k. In terms of these variables, the on-shell condition can 
be used to give 
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(B.16) 



with cos ipij = 1 — rii ■ rij, and </> the angle in A:j_-space (the azimuthal angle about Hi). We 
can do the k± and k + integrals using the on-shell and r a delta functions respectively. The 
resulting 5 , ™ eas (r*) has non-trivial integrals over k~ and cj): 
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Making the change of variables u = cot ^ -sjk+jkr = cot f ( J 2 °,we find that S^ eas (r^) 



can be written as 
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where I(a,/3,t) is defined in Eq. (B.10). Using Eq. (B.15) we find the result given in 
Eq. (5.18). 



B.2.3 SL 



The ©-functions in Eq. (5.14) are easiest to deal with if we shift to variables where each 
©-function is in a different variable. The simplest choices are just the arguments of the G 
functions Ga and @ l R , k° and u = cot ^\/k + jk~ , respectively, where k^ are defined with 
respect to direction rtj. This gives a form similar to the integral in Sfj eas (rl), 
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where I(ct,f3,t) is defined in Eq. (B.10) and evaluates to Eq. (B.15). This gives Eq. (5.15). 



B.3 ^ cas (r a fc ) and for fc 7^ t, j 

We again use light cone coordinates centered on jet k. The integrations involved in 
Sfj eas { T a) an d Sij only give rise to a 1/e pole as explained in the text, but integrating 
the eikonal factor l/(ra, • k)(rij ■ k) is more complicated than for the other cases since there 
is a third direction, ri/%, involved. 
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For unmeasured jets when there are n > 3 total final state jets, Sfj is needed. However, 
as we explain in the text, measured jets violate consistency at 0(l/t 2 ) even for n = 2 (non 
back-to-back) jets and the contribution of Sfj does not ameliorate this fact when n > 3. 
To show this, we need to evaluate the divergent contribution of Sfj. In addition, we give 
the form of the finite pieces which are 0(l/t 2 ). 

For each measured jet when there are n > 3, the sum 5'™ eas (r^) + Sfj8(r^) is needed. 
However, in this case the 1/e pole cancels in this sum and we are left with only a single, 
finite integral to evaluate. This is clear from the expressions for Sfj ea,s (T^) and which 
we derive in Sec. B.3.2 and Sec. B.3.3, respectively. We evaluate the sum explicitly in 
Sec. B.3.4. 



B.3.1 Common Integrals 

We find the following integral arising in both Sfj eas (r^) and S^: 
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where we defined 
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We can evaluate I^ ^ straightforwardly. For the range of our interest, t a>b > 1 and 
< u < 1, it gives 
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and 
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For r = 0, this simplifies to 
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Notice that both fi and / 2 are 0(l/t 2 ). 

The ©(e 1 ) piece, jW, is less trivial. However, the only property of 1^' that we need 
is that ^ 

Mt a ,t b ,P)= [ duuI ( > 1 \ U ;ta,t b ,(3) = 0(l/t 2 ), (B.28) 
Jo 

which can be seen by taking the large-i limit of I^ 1 ' in Eq. (B.21). The integral is finite an 
suppressed by 1/i 2 . 

B.3.2 S§ 

To compute S^, we choose a coordinate system such the is in the z-direction and Hi 
lies in the xz-plane. In terms of the light-cone coordinates about and the variable 
u = cot =j \Jk + jk~ , we have 
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where is defined as the angle between the ik- and j/c-planes. Using the relation 
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we find that can be written as 
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where I(u;t a ,h, f3ij) is defined in Eq. (B.20). Expanding in e, we find 
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B.3.3 5- cas (r a fc ) 

Using the same coordinate system as for S-?-, we find that 5*^ eas (r^) can be written as 
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and fx, /2, and / 3 are given in Eqs. (B.24), (B.25), and (B.28), respectively. 



B.3.4 S^ + S^(r k a ) 

The sum of Eqs. (B.32) and (B.35) is finite. We find 



a {T a ) + SijS{T a ) 



( 



4A 2 



tan' 



-2a 



i? 



+ (2 — a)f2(Uk,tjk, Pij, 0) — 2/2 (ii 2 , tife, /%, tan" —J 



2( ) fx(Uk,tjk, fiij) 



(B.37) 



where fx and / 2 are given in Eqs. (B.24) and (B.25), respectively. 



- 73 - 



C. Convolutions and Finite Terms in the Resummed Distribution 



In evaluating the final resummed distribution Eq. (6.40), each measured jet function must 
be convolved against a corresponding soft function piece 5' meas . These convolutions take 
the form. 



drjdrsdr'j dr' s J(t'j, tij)S™ as (T> s , Ms ) 



Q(rj ~ r'j) 
(rj - t'Y+^ 



e(r 5 - r' s ) 
( rs - T ')i+4 



5(t-tj-t s ) ■ 
(C.l) 

For the class of functions of the form x" 1 ^^ with ui ^ and ui < 1, we define the plus 
distribution by 



8(a) 



lim 

/3— >0 



@(x-/3) p 



-5{x- p) 



+£(- 



n=0 



Q(x) ln n x 



. (C2) 



where the plus functions on the second line are given by Eq. (A. 2), 

ln n+1 



~0(x)ln n (xy 


= lim 
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(C.3) 



From these definitions, we can derive the identities (see, e.g., Appendix B of [54]) 
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and 



, [ e(r-r') i , 
[{t-t') 1+uj \ + y 
Q(t-t') 1 re(r') 



r — r 



0(r - r' 



(r - r') 



r 

e^rQlnr 7 
r' 



e(r) 

Q(r) 



lnr-Hf-l-w] 



(C.5a) 
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Using the above identities, we find that the final result can be written in the form 



Eqs. (6.40) and (6.42) with the functions dj{r\) given by 
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and 
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for quarks and by 

d 9 f ™K) = C A 
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(C.6d) 



for gluons, where in all cases x con e, £1,2 (denned in Eqs. (A. 3) and (A. 8)) are evaluated at 
r = r j and r gj5 are denned in Eqs. (A. 13) and (A. 23). 
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D. Color Algebra for n = 2, 3 Jets 



For the two and three jet cases, there are no color correlations since all color generator 
inner products Tj • Tj can be expressed in terms of the Casimir invariants Ca and Cp- For 
n = 2, there is a quark jet with charge T q and an anti-quark jet with charge Tg that each 
square to Cf- There is only one inner-product in this case and using color conservation 
(Si = 0), we have that 

T g • T q = -T 2 q = -T? = -C F . (D.l) 
For n = 3 jets color conservation gives that, for example, 



1 



Ti • T 2 = t [( Tl + T 2 ) 2 - Tf - T|] 



= i[T|-Tf-T|] . (D.2) 
Referring to the quark, anti-quark, and gluon generators as T„, Tg-, and T s , respectively, 



using T 2 = T? = C> and T 2 = C A in Eq. (D.2) gives 

C A 



T g • T 5 = — 

T g • T g = Tg • T 9 = ■ (D.3) 
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